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1.  PROGRAM  OBJECTIVES 

The  goals  of  this  program  were  to  identify,  develop  and  implement  configurations  in  which  phonons  and  the  cou¬ 
pling  of  phonons  to  carriers  lead  to  enhanced  device  performance.  A  significant  component  of  our  efforts  centered 
on  artificial  semiconductor  heterostructures  such  as  quantum-wells  (QW),  -dots  (QD)  and  superlattices  (SL’s). 
Unlike  bulk  materials,  the  phonon  properties  of  these  low-dimensional  systems  and,  in  particular,  the  phonon  fre¬ 
quency,  group-velocity,  density  of  states  as  well  as  the  strength  of  the  interaction  with  carriers  can  be  widely  modi¬ 
fied  to  improve  the  performance  of  devices.  In  this  program,  we  considered  existing  devices  such  as  quantum- 
heterostructure  lasers  and  detectors  and  field-effect  transistors  (FET’s)  as  well  as  new  concepts  relying  on  the  inter¬ 
action  between  carriers  and  coherent  phonons,  specifically  phonon-assisted  lasing  and  tunneling.  Our  interactive 
theoretical  and  experimental  efforts  focused  on  (/)  the  development  of  models  to  predict  and  guide  the  experimental 
search  for  optimal  design  and  performance  of  various  electrical  and  optoelectronic  devices,  particularly  for  control 
of  thermal ization  rates,  (»)  the  improvement  of  existing  and  the  realization  of  new  coherent  optical  and  acoustic 
phonon  sources  using  both  ultrafast  lasers  and  electrical  methods,  and  ( iii )  the  application  and  improvement  of  state- 
of-the-art  materials  growth  and  micro-fabrication  methods  to  implement  new  ideas  for  control  of  phonon  parame¬ 
ters,  and  to  develop  phonon  filters,  lenses  and  reflectors. 

2.  PROGRAM  ACCOMPLISHMENTS 


2.1  Phonon  Generation  and  Detection 


2.1.1  COHERENT  PHONON-POLARITONS,  ACOUSTIC  AND  OPTICAL  PHONONS 


Phonon-Polaritons  -  Merlin  and  co-workers  used  an  optical  method  based  on  subpicosecond  pulses  moving  in  a 
nonlinear  medium  to  generate  Cherenkov  radiation  (CR)  in  the  infrared  range  through  phonon-assisted  phase  match¬ 
ing.  The  pulses,  of  central  frequency  <nL,  move  through  the  medium  with  velocity  v  =  cs(coL)  and  interact  with  them¬ 
selves  to  generate  a  low-frequency  polarization  through  frequency-difference  generation,  a  nonlinear  effect  associ¬ 
ated  with  the  susceptibility  tensor  x(2);  dajdq  and  q  are  the  light  group  velocity  and  wavevector.  The  generation 
of  CR  by  such  methods  was  pioneered  by  Auston  and  coworkers  [1]  who  applied  the  technique  to  LiTa03  in  a 
range  for  which  v  »  c0,  where  c0  is  the  speed  of  light  at  zero  frequency.  In  our  experiments  on  ZnTe  and  ZnSe,  we 
detected  emission  at  subluminal  speeds  involving  frequencies  below  those  of  the  corresponding  infrared-active  pho¬ 
non.  Results  for  for  ZnSe  were  reported  in  [2].  Data  for  ZnTe  is  shown  in  Fig.  1.  Compared  with  conventional  ul¬ 
trafast  methods  for  the  generation  of  coherent  THz  radiation,  our  approach  shows  some  advantages  on  how  to  attain 
phase-matching  conditions.  The  generation  of  Cherenkov 
radiation  by  electric  dipoles  generated  by  subpicosecond 
pulses  is  physically  indistinguishable  from  the  excitation 
of  coherent  polaritons  by  impulsive  stimulated  Raman 
scattering  (RS)  [2]. 


Acoustic  Modes  -  Bucksbaum  and  Merlin  [3]  proposed 
that  control  of  coherent  phonons  can  be  applied  to  slice  a 
beam  of  hard  x-rays  to  generate  subpicosecond  pulses. 
This  so-called  Bragg  switch  was  central  to  the  collabora¬ 
tive  effort  based  at  the  ANL  Advanced  Photon  Source  in¬ 
volving  MURI  researchers  and  colleagues  at  the  Univer¬ 
sity  of  Michigan  and  other  institutions.  These  efforts  led  to 
novel  ideas  for  the  development  of  ultrafast  x-ray  meth¬ 
ods,  as  well  as  real-time  studies  of  the  propagation  of  a 
high-frequency  sound  pulse  across  a  solid  [4].  Figure  2 
shows  the  time-resolved  transmission  for  sound  generated 
by  ultrafast  laser-excitation  at  t  =  0  on  the  exit  face  of  a  Ge 
crystal.  We  detected  both  the  forward  and  deflected  x- 
rays  as  a  function  of  time  and  crystal  angle.  The  x-rays 
propagate  via  the  Borrmann  effect.  Starting  at  t  =0  the  dif¬ 
fracted  intensity  switches  from  the  forward  beam  to  the 
deflected  beam  in  time  whose  measurement  is  limited  by 
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Fig.  1  -  Theoretical  and  experimental  plots  for  ZnTe. 
(a)  Calculated  E  vs.  p  and  (z/v-t)  for  a  point  dipole 
moving  at  v  ~  0.984c0.  (b)  Convolution  of  (a)  with  a 
Gaussian  function,  (c)  Differential  transmission  meas¬ 
urements  as  a  function  of  the  probe  beam  location.  For 
both  beams,  the  central  energy  and  the  width  of  the 
pulses  were  1.53  eV  and  -  50  fs. 
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the  x-ray  pulse  duration  (~  lOOps).  Ap¬ 
proximately  75%  of  the  energy  is  transferred 
from  one  beam  to  the  other.  Following  the 
initial  switch,  Pendellosung  oscillations  are 
apparent  in  both  the  forward  and  deflected 
beams.  The  oscillations  appear  to  be  n  out  of 
phase  with  a  frequency  of  ~1  GHz,  and  they 
decay  over  several  x-ray  attenuation  lengths. 
Around  60ns  the  acoustic  pulse  approaches 
the  input  face  of  the  crystal  at  the  sound 
speed,  so  that  the  acoustic  disturbance  forms 
a  thin  interface  separating  the  near-input  re¬ 
gion  from  the  rest  of  the  crystal.  This  has  a 
large  effect  on  the  coupling  of  the  input  beam 
into  the  a  mode  that  can  cross  the  thick  crys¬ 
tal  with  relatively  little  attenuation.  As  a  re¬ 
sult,  more  than  50%  of  the  input  intensity 
couples  to  the  a  wave. 

Superlattices  of  a  tailored  periodicity  provide 
a  means  to  couple  high-frequency  sound  to 
light  pulses.  As  opposed  to  picosecond  ultra¬ 
sonics  [5],  where  the  generated  acoustic 
wave  has  a  broad  frequency-spectrum  and  is 
usually  limited  to  a  few  hundreds  of  GHz, 
superlattices  can  act  as  sources  of  mono¬ 
chromatic  acoustic  waves  of  THz  frequen¬ 
cies.  Although  folded  acoustic  phonons  were 
discovered  nearly  twenty  years  ago,  the  has 
been  renewed  interest  due  to  the  observation 
of  coherently  generated  acoustic  phonons. 
Merlin’s  group  developed  a  model  for  under¬ 
standing  the  dynamics  of  the  high  frequency 
modes  in  layered  structures  and  superlattices. 
The  model  was  tested  by  performing  pump- 
probe  ultra-fast  optical  experiments  on  sam¬ 
ples  of  AlAs-GaAs  based  superlattices.  They 
also  developed  an  understanding  of  the  be¬ 
havior  of  the  modes  near  surfaces,  and  related 
these  modes  to  surface  avoiding  modes. 

Figure  3  illustrates  the  behavior  of  a  surface¬ 
avoiding  mode.  Calculations  of  the  normal 
modes  show  that  their  amplitude  is  a  maxi¬ 
mum  inside  the  superlattice  and  decreases 
towards  the  surfaces.  A  stationary  wave  can¬ 
not  satisfy  the  boundary  condition.  Since  the 
solutions  are  linearly  independent  at  the  zone 
center  (and  edge),  the  only  way  to  fulfill  the 
condition  is  to  avoid  the  surface  by  making 
the  amplitude  of  the  wave  close  to  zero  in  its 
vicinity.  Such  modes  are  known  as  surface 
avoiding  modes. 
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Fig.  2  -  Time-resolved  diffraction  reflection  following  laser  impulsive 
excitation  of  an  acoustic  wave  on  the  exit  face  of  a  thick  Ge  [001] 
crystal.  Top:  Forward  (blue)  and  deflected  (red)  x-ray  intensity  as  a 
function  of  time  following  excitation.  Bottom:  expanded  scale  show¬ 
ing  the  regions  where  the  acoustic  pulse  is  near  the  exit  or  entrance 
crystal  surface.  2-d  plots:  x-ray  transmission  vs.  time  and  crystal  an¬ 
gle  near  the  forward  (top)  and  deflected  (bottom)  diffraction  peaks. 


z  [nm] 


Fig.  3  -  Calculated  dispersion  of  LA  modes  propagating  along  [001]  for 
the  GaAs-AlAs  SL  investigated  in  this  work.  Panels  on  the  right  show  a 
schematic  diagram  and  four  representative  eigenmodes  of  the  619-nm- 
thick  SL  slab  (striped  area)  attached  to  a  5-pm-thick  substrate.  Vertical 
green  bars  denote  the  SL-substrate  boundary.  SAM  is  the  surface¬ 
avoiding  mode  at  q  ~  0.  The  wavevector  associated  with  modes  labeled 
FP  is  sufficiently  far  removed  from  the  zone  center  so  that  the  FP-waves 
reflect  only  weakly  at  the  boundary.  The  mode  with  frequency  in  the 
minigap,  GM,  decays  exponentially  into  the  SL. 


Time-domain  experiments  were  performed  on  a  superlattice,  grown  by  molecular  beam  epitaxy,  of  75  periods  of  59 
A  GaAs  and  23.5  A  AlAs.  A  70  nm  layer  of  AlAs  was  grown  between  the  superlattice  and  the  substrate  to  act  as  an 
etch  stop.  The  time-resolved  experiments  were  carried  out  in  the  standard  pump  and  probe  configuration,  using  an 
optical  parametric  amplifier  which  delivered  ultrafast  pulses  (70  fs)  at  a  rate  of  250  KHz  tunable  in  the  visible  range 
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(400-650  nm)  with  a  pulse  energy  of  1  nJ.  These  pulses  '  1  ’  1  ’  1  ’  1  ’  1  ’ 

are  split  into  pump  and  probe.  The  pump  excites  co-  I.OxlO'5  -  I 

herent  vibrations  of  the  lattice  that  modulate  the  index  f| 

of  refraction  and,  hence,  the  reflectivity  of  a  time-  5.0x10  -  r  ,i  ii  i  '•••.  ,  i 

delayed  probe  pulse.  The  pulses’  central  wavelength  is  q-  q q  ,  Ml  jra  t  i  JL  i  uJyi  ■ 

tuned  for  resonant  excitation  with  the  superlattice  at  &  if  IJ  |i  If  \/¥  '  *  "  "x,  ^W. 

530  nm.  The  pump  beam  is  chopped  at  2  KHz  and  the  -5.0x1  O'*  -  II  f  »  |  I  |  v 

signal  is  analyzed  with  a  lock-in  amplifier  at  this  fre-  R  _ pj — - — — - . — Xi  ■ 

quency.  All  the  experiments  were  done  at  room  tern-  -1.0x10  -  \  (uL  , 

perature.  In  Fig.  4  we  show  the  differential  reflectivity  -1 5xio'5  -  \ 

data  as  a  function  of  the  probe  time  delay  between  0 

and  250  ps  after  subtracting  a  slowly  varying  back-  -2.0x1 0'5  — — i - >■  — ;  -----  ,  — — — 

ground.  Similar  scans  were  taken  up  to  500  ps  time  de-  0  50  100  150  200  250 

lay.  An  oscillation  with  a  period  of  14  ps  can  be  seen  time  delay  [ps] 

and,  superimposed,  higher  frequency  oscillations  as 

shown  in  the  inset.  Note  that  what  may  appear  to  be  Fig.  4  -  Differential  reflectivity  as  a  function  of  time  delay, 
noise  is  actually  the  superposition  of  modes  of  differ¬ 
ent  frequencies.  The  Fourier  transform  shows  clear  peaks  even  after  time  delays  of  800  ps.  The  fast  Fourier  trans¬ 
form  of  the  time  trace  in  Fig.  4  is  shown  in  Fig.  5  together  with  the  dispersion  relation  shown  in  Fig.  3.  Four  peaks 
can  clearly  be  seen  in  the  power  spectrum.  By  studying  the  Fourier  transform  of  different  time  sections  we  can  ob¬ 
tain  information  on  the  evolution  of  the  four  different  modes.  In  Fig.  5  we  show  the  Fourier  transform  of  time  scans 
similar  to  the  one  shown  in  Fig.  4  for  four  different  time  intervals.  At  short  times,  four  peaks  are  seen.  Here,  the  trip¬ 
let  near  0.6  THz  consists  of  two  peaks  corresponding  to  the  modes  for  q  =  2n A,  and  a  peak  slightly  shifted  from  the 
center  of  the  other  two  that  corresponds  to  one  of  the  modes  at  the  center  of  the  mini-Brillouin  zone.  As  we  move  to 
long  times,  spectra  (b)  and  (c),  the  doublet  (peaks  2  and  4)  fades  away  and  only  the  mode  at  q  =  0  remains  visible  in 
the  spectrum  together  with  the  Brillouin  peak.  At  later 
times,  spectrum  (d),  we  notice  that  the  zone  center 

mode  is  still  visible  but  the  Brillouin  peak  has  almost  — — r — , — , — , — , — , — , — , — , — , — , — T _ , — , — r — 

disappeared.  The  confined  nature  of  this  mode  mani¬ 
fests  itself  in  time-resolved  experiments  as  an  oscilla-  i  ■  A  - 

tion  with  a  very  long  lifetime.  Calculations  of  the  ^  - — 

acoustic  normal  modes  show  that  the  mode  with  the  -  i 

very  long  lifetime  is  indeed  localized  in  the  superlattice.  §  I 

Finally  we  note  that  the  apparent  short  lifetime  of  the  £*  i  — • - — r  \ — 

traveling  modes  is  due  to  their  disappearance  from  the  “  j 

surface,  rather  than  decay  of  the  phonon  into  other  f  y  ,  i  ^ 

modes.  <u  K  V \  aA\  A'-  a- vaJ  yv.A 

Optical  Phonons  -  While  there  has  been  widespread  —  !  \  3 

agreement  that  stimulated  RS  is  the  key  generation  •  i  "I  2  1  4 

mechanism  for  coherent  phonons  in  the  transparent  re-  ]  \i\i.  ftil  a  a 

gion,  its  role  in  opaque  materials  has  remained  contro-  -  j  ^  -  « 

versial.  Merlin  and  co-workers  demonstrated  a  previ-  '‘j’  ‘  ‘  j  '  '  \  |  1  \  \  |  j-1-!  . : 

ously  unrecognized  feature  of  RS  that  settles  the  discus-  3  /  \  / 

sion  [6].  Specifically,  they  proved  that  stimulated  RS  is  ^  j  /  \  / 

defined  not  by  one,  but  by  two  distinct  tensors.  The  first  7)  2  \  /  \  j  / 

one  is  the  standard  Raman  susceptibility  discussed  in  |  /  \  / 

textbooks  which  gives  both  the  modulation  of  the  linear  |  i  -  j/ _ \  i  / _ ; 

susceptibility  by  the  phonon  and  the  cross  section  for  %  /  V  / 

spontaneous  RS.  The  second,  new  tensor,  characterizes  oo  "  0I2 . 0.3  '  0I4  '  0)5  0.6  '  oj  '  ole  '  09 

the  electrostrictive  force  acting  on  the  ions.  Consistent  Frequency  [THz] 

with  the  common  understanding  of  stimulated  Raman 

processes,  the  real  parts  of  these  tensors  are  identical  Fig.  5  -  Fourier  transform  of  different  time  sections  like  the  one 
and,  thus,  they  become  one  and  the  same  in  the  domain  shown  in  Fig.  4  together  with  the  dispersion  relation  calculated 
of  transparency.  Relevant  to  the  absorbing  region,  how-  with  the  parameters  of  the  sample.  The  spectra  correspond  to, 
ever,  their  imaginary  components  differ  appreciably,  a  =  0-125  ps,  b=  125-250  ps,  c  =  250-350  ps  and  d  =  350-500  ps. 


Frequency  [THz] 

Fig.  5  -  Fourier  transform  of  different  time  sections  like  the  one 
shown  in  Fig.  4  together  with  the  dispersion  relation  calculated 
with  the  parameters  of  the  sample.  The  spectra  correspond  to, 
a  =  0-125  ps,  b  =  125-250  ps,  c  =  250-350  ps  and  d  =  350-500  ps. 
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These  predictions  were  tested  by  examining  the  behavior  of  coherent  phonons  generated  at  photon  energies  in  the 
range  of  the  E2  critical  point  of  Sb.  The  experimental  results  are  in  very  good  agreement  with  the  two-tensor  model. 

The  anharmonic  properties  of  the  low-frequency  E2  phonon  in  ZnO  were  measured  by  Merlin  and  co-workers  using 
impulsive  stimulated  RS  [7].  At  5  K,  the  frequency  and  lifetime  are  (2.9789  ±  0.0002)  THz  and  (211  ±  7)  ps.  The 
unusually  long  lifetime  and  the  high  accuracy  in  the  determination  of  the  frequency  hold  promise  for  applications  in 
metrology  and  materials  characterization.  The  temperature  dependence  of  the  lifetime  is  determined  by  two-phonon 
up-conversion  decay  contributions,  which  vanish  at  zero  temperature.  Results  shown  in  Fig.  6  suggest  that  the  life¬ 
time  is  limited  by  isotopic  disorder  and  that  values  in  the  nanosecond  range  may  be  achievable  in  isotopically-pure 
samples.  Everitt’s  group  also  investigated  coherent  optical  phonons  in  bulk  ZnO  using  samples  that  were  thermally 
treated  under  forming  gas  environment  (from  Cermet,  Inc.  and  Eagle-Picher).  Phonon  lifetimes  can  be  extracted 
from  impulsive  stimulated  RS  measurements.  As  seen  in  Fig.  7,  the  preliminary  measurements  at  Duke  University 


lineshape  FWHM  from  Lorentzian  fits  to  the  fre¬ 
quency-domain  data  (circles)  and  time-domain  fits  to 
the  coherent-oscillation  amplitude  using  a  simple  de¬ 
caying  exponential  (squares).  Full  and  dotted  lines  are 
fits  neglecting,  respectively,  down-  and  up-conversion 
processes.  Inset:  Frequency  vs.  temperature.  The  dot¬ 
ted  curve  is  a  guide  to  the  eye. 


time  delay  (ps) 

Fig.  7  -  Impulsive  stimulated  Raman  data  for  a  bulk  ZnO  sample 
(Cermet).  The  inset  shows  the  fast  Fourier  transform  of  the  data. 


using  a  mode  locked  Ti:Sapphire  laser  at  800  nm  with  ~100  fs  pulses  show  a  reasonable  ISRS  signal  at  room  tem¬ 
perature.  The  inset  shows  the  fast  Fourier  transform  of  the  data,  which  has  two  well  defined  peaks  at  -2.05  and 
-2.75  THz.  Lee  et  al.  [8]  identified  the  Ai(LO),  E2-high,  and  E2-I0W  optical  phonon  mode  frequencies  as  17.2, 
13.2,  and  2.98  THz,  respectively,  for  a  hydrothermal ly-grown  bulk  ZnO  sample.  For  the  observation  of  E2-high,  and 
E2-I0W  modes,  pump  and  the  probe  beams  need  to  be  orthogonal;  however,  the  A](LO)  mode  may  be  observed  in 
any  geometry.  In  these  measurements,  the  pulse-width  of  -  lOOfs  limited  the  observation  of  any  modes  with  fre¬ 
quencies  above  1 0  THz.  The  observed  frequency  of  2.75  THz  is  very  close  to  the  value  reported  for  the  E2-low 
mode  by  Lee  et  al.  The  dephasing  time  at  room  temperature  is  less  than  10  ps. 

Acoustic  Modes  in  Quantum  Dots  -  Studies  of  acoustic  modes  in  nanoparticles  were  pursued  in  Merlin’group  [9]. 
Using  ultrafast  optical  pulses,  coherent  phonons  were  generated  in  a  semiconductor  doped  glass,  made  of  semicon¬ 
ductor  inclusions  of  nanometer  size  in  an  insulating  matrix.  By  tuning  the  laser  to  resonate  with  the  absorption  edge 
of  the  semiconductor  QDs,  they  were  able  to  generate  high  amplitude  coherent  optical  and  acoustic  phonons.  The 
stimulated  RS  model  reproduces  accurately  the  experimental  results  as  to  be  the  mechanism  responsible  for  the  pho¬ 
non  generation  in  the  absorbing  regime.  These  studies  showed  that  the  initial  phase  of  the  coherent  field  carries  in¬ 
formation  about  coupling  between  modes.  If  the  modes  are  coupled,  the  initial  phase  reflects  this  effect  by  a  phase- 
shift  respect  to  the  phase  predicted  for  independent  oscillators.  They  experimentally  detected  this  phase-shift,  due  to 
the  coupling  between  two  optical  modes,  in  the  time-domain  measurements.  At  low  power  excitation,  they  observed 
that  the  frequency  of  the  optical  coherent  phonons  deviates  from  values  obtained  from  spontaneous  Raman  scatter¬ 
ing,  as  shown  in  Fig.  8.  This  effect  was  ascribed  to  the  presence  of  electronic  impurity  states  (traps)  which  modify 
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Power  density  (W/cmJ) 

Fig.  8  -  Power  dependence  of  the  frequency  of  the  CdSe-like 
LO  phonon  for  Raman  and  pump-probe  measurements. 


Fig.  9  -  (left)  Energy  diagram  of  the  QD  with  a  single  defect, 
(right)  Time  dependence  of  the  population  of  the  lowest  level. 


the  nanocrystal  dielectric  function  in  the  infrared  and,  thereby,  renormalizes  the  frequency  of  the  optical  modes;  see 
the  energy  level  diagram  in  Fig.  9(left).  In  particular,  a  filled  electronic  trap  with  an  allowed  transition  frequency 
larger  than  the  bare  optical  phonon  frequency,  will  push  the  effective  optical  phonon  frequency  detected  to  lower 
values.  This  happens  for  all  the  power  densities  used  in  the  cw  experiment,  for  which  the  trap  is  always  full.  In  time- 
domain  experiments,  this  effect  operates  only  if  the  trap  filing  time  is  smaller  than  the  phonon  time  decay,  as  t\<rph 
in  Figure  9  (right),  which  correspond  to  power  P,  (n0J  is  the  population  of  the  trap  level  0d).  If  the  power  density  is 
always  higher  than  the  power  to  saturate  the  trap,  P0,  it  was  shown  that  moving  the  power  density  in  time  domain 
experiments,  it  is  possible  to  go  from  the  renormalized  frequency  (at  power  /“,)  to  the  bare  frequency  (at  power  P2), 
which  is  consistent  with  the  behavior  seen  in  Fig.  8. 


2.1.2  ACOUSTIC  PHONONS  IN  THE  NITRIDE  SYSTEM 

The  immense  potential  of  acousto-optical  devices  based  on  nitrides  as  well  as  the  poorly  understood  physics  of 
acoustic  phonons  in  multiple  quantum  wells  (MQWs)  make  the  study  of  acoustic  phonons  particularly  important. 
The  collaboration  between  the  groups  of  Hadis  Morkof  at  Virginia  Commonwealth  University,  and  Henry  Everitt,  at 
Duke  University,  focused  on  the  generation  of  coherent  acoustic  phonons  in  nitride  superlattices.  Using  sub¬ 
picosecond  optical  pump-probe  techniques,  Everitt  and  coworkers  were  able  to  generate  and  control  coherent  zone- 
folded  longitudinal  acoustic  phonons  (ZFLAPs)  in  InGaN  multiple  QW  structures  [10],  Using  differential  transmis¬ 
sion  (DT)  techniques,  their  measurements  revealed  that  carriers  injected  near  the  barrier  band  edge  were  quickly 
captured  into  the  quantum  wells  and  generated  strong  coherent  ZFLAP  oscillations  of  frequency  0.69  THz.  Two- 
pump  DT  was  used  to  generate  coherent  ZFLAP  oscillations  whose  relative  timing  and  amplitude  could  be  con¬ 
trolled.  Enhancement  and  suppression  of  ZFLAP  oscillations  were  demonstrated.  The  ability  to  generate  and  control 
acoustic  phonons,  dramatically  simplified  by  the  properties  of  the  InGaN  MQW  system,  is  of  scientific  and  techno¬ 
logical  utility.  By  canceling  coherent  ZFLAPs  one  half  period  after  their  creation,  a  single  mode  acoustic  phonon 
impulse  is  generated  which  can  propagate  through  the  sample  and  whose  subsequent  effects  on  the  material  response 
can  be  monitored.  It  should  be  possible  to  probe  the  effects  of  single  mode  acoustic  phonon  interactions  with  carriers 
or  excitons  in  a  controlled  manner.  More  generally,  complex  temporal  pump  pulse  waveforms,  created  by  genetic 
learning  algorithms  and  spatial  light  modulators,  can  be  used  to  create  non-trivial  phonon  excitations  to  explore  or 
control  phonon-mediated  decay  or  dephasing.  In  addition,  spectrally  pure,  coherent  0.10  THz  bulk  and  0.69  THz 
zone  folded  acoustic  phonons  were  generated  using  similar  impulsive  techniques.  As  before,  a  single,  high  fre¬ 
quency  acoustic  phonon  mode  was  excited  using  impulsive  differential  transmission,  and  it  was  found  that  these 
zone  folded  acoustic  phonons  are  insensitive  to  temperature  changes  between  80K  and  300K.  Acoustic  modes  at 
both  frequencies  were  simultaneously  excited  using  impulsive  differential  reflection.  It  was  observed  that  the  bulk 
phonons  travel  ballistically  through  the  cap  layer  and  reflect  efficiently  from  the  semiconductor/air  interface, 
whereas  reflected  zone-folded  phonons  were  not  observed. 

Since  the  coherent  acoustic  phonon  generation  is  directly  linked  to  the  carrier  dynamics,  it  was  important  to  under¬ 
stand  the  relaxation  pathways  including  carrier  capture  and  recombination  Recent  measurements  of  ultrafast  dynam- 
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ics  in  InGaN  and  AlGaN/GaN  MQWs  helped  to  identify  different  processes  from  sub-picosecond  to  nanosecond 
timescales  [11].  Carrier  capture  into  the  QWs  is  observed  to  occur  in  less  than  1  ps  and  this  resulted  in  an  immediate 
relief  of  the  strain  used  to  generate  coherent  acoustic  phonons. 

Another  important  objective  of  the  Morkof-Everitt  collaboration  was  to  build  an  acoustic  phonon  resonator  with  a 
frequency  of  ~  1  THz.  To  make  this  achievable,  the  coherence  time  of  the  acoustic  phonons  needed  to  be  increased 
and  the  leakage  from  the  MQW  region  had  to  be  mini¬ 
mized.  For  this  reason,  AlN/GaN  MQW  layers  were 
grown  by  plasma  assisted  MBE  to  be  used  as  phonon 
mirrors  on  both  sides  of  an  InGaN  MQW  structure.  The 
QW  periods  were  determined  from  high-resolution  x- 
ray  measurements.  The  superlattice  satellite  peaks  indi¬ 
cate  a  periodic  structure  with  smooth  interfaces. 

For  acoustic  phonon  generation,  InGaN  MQW  struc¬ 
tures  were  grown  by  MOCVD  on  c-plane  sapphire  sub¬ 
strates  with  ~3  pm  GaN:Si  (1018  cm'3  doping)  template 
layers.  Three  samples  (CVD58,  CVD59,  CVD75)  with 
10  periods  of  In003GaN/In0  isGaN  were  characterized 
by  x-ray  diffraction,  photoluminescence  (PL),  PL  exci¬ 
tation  (PLE),  and  absorption  measurements.  One  of  the 
samples  (CVD75)  was  capped  with  a  ~50  nm  GaN 
layer.  PL,  PLE,  and  absorption  data  are  shown  in 
Fig.  11.  All  three  samples  exhibited  intense  lumines¬ 
cence.  However,  the  barrier  absorption,  which  was 
shown  previously  to  be  effective  in  acoustic  phonon 
generation,  is  low.  This  explains  why  time-resolved 
pump/probe  (differential  transmission)  measurements 
did  not  show  any  sign  of  acoustic  phonon  generation. 

Research  is  still  underway  to  optimize  the  growth  con¬ 
ditions  of  the  InGaN  MQWs  to  increase  the  barrier  ab¬ 
sorption  and  to  obtain  better  interfaces.  These  MQW 
structures  will  be  grown  on  GaN  templates  with  SiN 
nanonet  and/or  AlN/GaN  superlattice  structures  to  re¬ 
duce  the  defect  density  in  the  MQW  region.  After  the 
InGaN  MQW  quality  is  optimized,  same  structures  will 
be  grown  on  AlN/GaN  MQW  phonon  mirrors. 

2.1.3  ELECTRICAL  GENERATION  AND  AMPLIFICATION  OF  COHERENT  PHONONS 

Phonon  Generation  in  Nanoscale  Structures  Under  the  Overshoot  Transport  Regime  —  Kim’s  group  analyzed  the 
interaction  of  high-speed  electrons  and  polar  optical  phonons  in  nanoscale  semiconductor  structures.  Under  the  ve¬ 
locity  overshoot  regime,  it  is  found  that  electrons  can  induce  optical  phonon  instability;  i.e.,  the  lattice  optical  vibra¬ 
tions  grow  in  time  when  the  electric  current  exceeds  a  threshold  value.  For  nanoscale  diodes  based  on  III-V  com¬ 
pounds,  the  threshold  current  density  is  estimated  to  be  of  the  order  of  hundred  kA/cm2  and  the  corresponding  bias 
0.5-1. 5  V  for  the  diode  lengths  of  50-100  nm;  see  Fig.  12.  The  wavelengths  of  unstable  phonons  are  approximately 
10  nm  with  the  characteristic  time  for  instability  development  in  the  10  ps  range. 

THz  Microwave  Emission  in  Ballistic  Structures  Under  Optical  Phonon  Generation  -  In  nanostructures  with  the 
ballistic  electron  transport,  we  found  the  frequency  windows  with  negative  dynamic  resistance  due  to  the  effect  of 
generation  of  optical  phonons.  These  windows  are  near  the  optical  phonon  frequency.  Kim’s  group  calculated  the 
impedance  of  the  ballistic  devices  and  established  that  the  device  can  be  used  for  generation  of  THz  emission.  The 
typical  parameters  obtained  for  the  devices  based  on  III-V  compounds  are:  the  device  length  of  about  50-70  nm,  the 
electric  bias  of  0. 1-0.2  V,  and  the  frequencies  of  7-8  THz.  The  low  temperature  is  assumed  for  ballistic  transport. 

Generation  of  High-Frequency  Acoustic  Phonons  by  ac  Fields  in  Acoustic  Resonators  -  Kim’s  group  proposed 
and  analyzed  generation  of  THz  acoustic  phonon  pulses  by  using  the  ac  (MHz  range)  electric  field.  The  device  under 
consideration  is  a  n+(3D)-n(2D)-n+(3D)  or  a  p+-p-p+  diode.  The  base  of  the  diode  consists  of  a  heterostructure  with 
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Fig.  10  -  Proposed  sample  structure  to  increase  the  confine¬ 
ment  of  coherent  acoustic  phonons  in  the  InGaN  MQW. 
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Fig.  1 1-  PL,  PLE,  and  absorption  (ABS)  data  for  three 
InGaN  MQW  samples. 
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a  single  or  multiple  quantum  wells,  which  confine 
both  electrons  and  phonons.  The  lateral  transport  in 
the  base  with  2DEG  provides  amplification  of  a  pulse 
(packet)  of  the  THz  phonons.  The  pulse  moves  back 
and  forth  reflected  at  the  surfaces  of  the  diode.  Gen¬ 
eration  of  the  THz  phonons  arises,  when  the  ac  elec¬ 
tric  field  is  synchronized  with  the  phonon  pulse  mo¬ 
tion.  For  example,  the  generator  utilizing  a  p-doped 
SiGe  heterostructure  with  the  base  length  of  50.100 
pm,  can  generate  pulses  of  acoustic  phonons  in  the 
0.3-0.4  THz  range  by  using  the  electric  fields  of  0.5-1 
kV/cm  with  the  ac  frequency  of  20-50  MHz. 

Light  Modulation  by  High-Frequency  Phonons  - 
Using  a  photoelasticity  model,  we  studied  light 
modulation  by  coherent  optical  phonons  generated  in 
AlGaAs/GaAs  superlattices.  For  the  case  of  a  plane 
wave  case,  the  amplitude  of  the  scattered  light  wave 
was  estimated  to  be  approx.  HF4  to  10~3  of  the  inci¬ 
dent  amplitude.  A  further  analysis  demonstrated  that 
the  efficiency  of  light  modulation  can  be  enhanced 
drastically  if  the  phonon  generator  is  placed  inside  a 
light  cavity.  In  this  case  of  high-Q  Fabry-Perot  cav¬ 
ity,  the  amplitude  of  the  scattered  wave  can  be  com¬ 
parable  to  that  of  the  incident  wave.  This  suggests  that  the  phonon  generators  combined  with  the  optical  cavity  can 
be  used  for  very  efficient  cw  THz  modulation  of  light,  which  is  not  available  at  present  by  other  means. 

Phonon  Generation  by  Carrier  Drift  and  Cherenkov  Emission-  Kim  and  collaborators  developed  a  model  for  the 
amplification  of  confined  vibrational  modes  by  the  drift  of  the  two-dimensional  carriers  as  a  function  of  phonon  fre¬ 
quency,  temperature  and  heterostructure  parameters.  Two  electron-phonon  interaction  mechanisms  were  taken  into 
account  involving  interaction  via  the  deformation  potential  (p-SiGe/Si  and  n-AIGaAs/GaAs)  and  the  piezoelectric 
interaction  (n-AIGaAs/GaAs).  It  was  found  that  an  amplification  coefficient  of  the  order  of  hundreds  of  cm'1  for  the 
AlGaAs  heterostructures  and  thousands  of  cm"1  for  the  SiGe  heterostructures  can  be  obtained  in  spectrally-separated 
narrow  amplification  bands.  Amplification  of  shear  vertical  waves  in  the  heterostructures  grown  in  the  direction  of 
high  symmetry  (for  instance,  [001])  is  found  to  be  significantly  stronger  than  that  predicted  by  the  model  of  elasti¬ 
cally  isotropic  media.  The  NCSU  team  used  the  Cherenkov  generation  model  to  investigate  optical  phonon  amplifi¬ 
cation  by  a  drifting  degenerate  electron  gas.  Extension  of  the  analysis  for  the  degenerate  gas  is  important  since  this 
is  the  case  most  favorable  to  meet  the  generation  criteria.  The  results  demonstrate  high  efficiency  of  optical  phonon 
generation  as  well  as  its  significant  impact  on  the  electron  transport  in  the  generation  regime.  The  possibility  of  op¬ 
tical  phonon  avalanche  was  also  investigated.  A  comparison  between  the  experimental  data  (provided  by  Prof.  K. 
Tsen  of  Arizona  State  University)  and  the  theoretical  calulations  shows  a  good  fit  across  the  wavevector  range  under 
consideration  including  the  presence  and  location  of  the  threshold  wavevector.  This  observation  clearly  indicates 
that  the  unusally  large  optical  phonon  population  obtained  in  the  light  scattering  measurements  is  due  indeed  to 
Cherenkov  amplification  of  nonequilibrium  phonons. 

2.1.4  IMAGING  NANOSTRUCTURES  WITH  ACOUSTIC  PULSES 

Norris’  group  developed  a  novel  acoustic  imaging  scheme  based  on  ultrafast  optical  excitation  and  detection  of  high 
frequency  ultrasound  (coherent  longitudinal  phonon  pulses).  The  basic  idea  of  the  experiment  is  illustrated  in 
Fig.  13.  The  sample  is  a  0.5  mm  Si  wafer  coated  on  both  sides  with  thin  A1  films,  one  of  which  is  lithographically 
patterned  with  1-micron-width  lines  and  grids.  The  pump  pulse  heats  the  surface  of  the  film  which  causes  a  thermal 
expansion  which  in  turn  launches  an  acoustic  strain  pulse  (coherent  longitudinal  phonons)  into  the  film  and  sub¬ 
strate.  The  single  cycle  pulse  so  generated  will  have  a  peak  frequency  and  bandwidth  both  ~  100  GHz  and  will  be  a 
few  picoseconds  in  duration.  The  pulse  propagates  across  the  0.5  mm  Si  wafer  at  a  speed  vL,  equal  to  the  velocity  of 
a  longitudinal  sound  wave  in  Si  (1 1 1),  until  it  reaches  the  A1  film  on  the  opposite  side.  The  change  in  reflectivity 
caused  by  the  strain  pulse  can  be  detected  there  by  a  time  delayed  optical  probe  pulse  derived  from  the  same 
Ti:sapphire  laser  as  the  pump.  Since  the  characteristic  diffraction  length  for  a  single  cycle  pulse  is  given  by  the  ex- 


Fig.  12  -  Current  induced  increment  F„  for  different  optical 
modes  in  a  75-nm  long  GaAs  n-i-n  diode  in  the  overshoot  trans¬ 
port  regime.  The  current  7=1  corresponds  to  33  kA/cm2.  For 
7=  20,  the  applied  voltage  is  0.55  V.  The  inset  shows  Tn  vs  the 
mode  number  n  at  7  =  4.  The  dashed  line  indicates  the  optical 
phonon  loss  at  T=  100  K. 
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Fig.  13  -  Schematic  diagram  of  the  experiment. 


pression  z  =  r*/vt0 ,  where  v  is  the  velocity  (9355  m/s)  and  t0  is 
the  temporal  pulse  width  (~2.5  ps)  they  find  that  for  a  radius, 
r0  ~  3pm ,  z  ~  0.4  mm  which  is  less  than  the  thickness  of  the 
wafer.  Therefore,  it  is  possible  to  probe  the  picosecond  acous¬ 
tic  pulse  in  the  far  field  regime  for  a  tightly  focused  pump  spot, 
and  the  near  field  regime  for  a  larger  pump  spot. 

It  is  important  to  note  that  in  this  proof-of-concept  experiment, 
the  object  to  be  imaged  is  actually  the  source  of  the  detected 
acoustic  waves.  The  purpose  of  such  a  relatively  thick  sample  is 
to  allow  the  scattered  acoustic  pulses  to  propagate  a  sufficient 
distance  so  that  variations  in  the  spatial  profile  of  the  measured 
acoustic  field  can  be  probed  with  a  diffraction-limited  optical 
probe  spot.  This  long  acoustic  propagation  length  requires  the  sample  to  be  held  at  cryogenic  temperatures  in  order 
to  suppress  the  phonon-phonon  scattering  that  strongly  attenuates  acoustic  waves  of  such  high  frequency.  The 
changes  in  reflectivity  caused  by  the  propagating  phonon  pulses  are  then  measured  with  time-delayed  optical  probe 
pulses,  giving  measurements  of  the  scattered  acoustic  field.  The  measured  field  is  then  used  as  the  input  to  a  “time- 
reversed”  Huygens-Fresnel  diffraction  equation  in  a  numerical  backpropagation  algorithm  in  order  to  generate  a  re¬ 
construction  of  the  Al  pattern  that  launched  the  waves. 

Figs.  14  (a)  and  14  (b)  show  the  reconstructed  acoustic  image  of  an  array  of  1  micron  lines  and  an  SEM  of  the  same 
structure,  respectively  (scale  bar  on  all  figures  =  5  pm).  To  obtain  this  image,  it  was  necessary  to  account  for  the 
elastic  anisotropy  of  the  crystalline  Si  wafer.  The  image  shown  in  Fig.  14  (c)  is  not  corrected  for  this  anisotropy  ef¬ 
fect.  It  is  evident  from  the  figures  that  the  periodicity  and  shape  of  the  acoustic  image  was  incorrect  prior  to  the  use 
of  an  anisotropic  sound  velocity  calculated  from  the  known  elastic  constants  of  Si. 

After  establishing  the  correct  backpropagation  algorithm  for  the  acoustic  imaging  scheme,  the  first  2-D  ‘picosecond 
ultrasonic’  images  of  nanostructures  were  obtained  [12].  Fig.  14  (d)  shows  a  lineout  from  the  image  of  Fig.  14  (a), 
and  illustrates  that  the  imaging  resolution  in  these  figures  is  a  modest  700  nm,  far  from  the  optimum  -30  nm  resolu¬ 
tion  that  could  theoretically  be  achieved  with  such  high  frequency  acoustic  waves.  Future  plans  include  resolution 
improvement  efforts,  as  well  as  making  the  first  images 
of  subsurface  structures,  which  will  showcase  the  true 
advantage  of  an  ultrasonic  technique.  Related  work  will 
include  the  investigation  of  different  materials  as  trans¬ 
ducer  layers,  in  the  hopes  of  improving  fundamental  sig- 
nal-to-noise  level  in  the  detection  of  small  acoustic  tran¬ 
sients. 

In  addition  to  imaging  with  coherent  phonons,  Norris’ 
group  investigated  phonon  pulse  propagation  (disper¬ 
sion,  diffraction,  and  nonlinearity  all  must  be  consid¬ 
ered)  with  a  similar  experimental  arrangement.  This 
work  is  of  interest  both  from  a  fundamental  standpoint 
as  well  as  through  its  usefulness  in  determining  the  vi¬ 
ability  of  imaging  schemes.  The  key  contribution  made 
through  this  investigation  was  to  quantify  the  effects  of 
diffraction  on  coherent  phonon  pulses  in  an  anisotropic 
crystalline  medium.  Their  work  in  this  direction  led  to  a 
fruitful  collaboration  with  the  theoretical  group  of  Jacob 
Khurgin  at  Johns  Hopkins  University. 


Fig.  14  -  Phonon  imaging  of  1  pm  lines. 


2.1.5  PHONONIC  IMAGING  OF  INTERFACES 

Asymmetric  interfaces  exhibit  both  mode-conserving  and  mode-converting  scattering.  Thus,  they  display  the  full  po¬ 
tential  of  tri-refringent  acoustic  materials.  This  complexity  could  be  an  advantage  for  phonon  optics  applications. 
For  example,  a  well  designed  asymmetric  interface  could  be  used  to  convert  the  longitudinal  waves  normally 
launched  by  picosecond  excitation  of  heterostructures  to  transverse  bulk  waves.  The  ability  to  generate  both  longitu¬ 
dinal  and  transverse  bulk  waves  would  greatly  extend  the  current  range  of  picosecond  ultrasonic  experiments. 
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Msall’s  work  on  a  grain  boundary  in  CaF2  highlights  the  challenges  of  creating  detailed  models  of  phonon  scattering 
at  asymmetric  interfaces.  High  quality  modeling  of  asymmetric  interfaces  requires  extremely  accurate  structural  in¬ 
formation  available  that  we  obtain  using  Electron  Backscatter  Difffactometry  (EBSD).  Fig.  15  shows  the  phonon 
flux  pattern  for  a  CaF2  sample  in  which  the  first  section  is  characterized  by  Euler  angles  91=  2°,  0i  =  17°,  91  =  37°, 
and  the  second  by  92=  57°,  02  =  36°,  v|/2  =  62°.  The  features  in  the  computer  simulation  (Fig.  15  (a)  and  (b))  compare 
well  with  those  in  the  experimental  image,  Fig.  15  (c).  There  are  minor  distortions  in  some  feature  shapes  (e.g.,  the 
feature  marked  5)  because  the  grain  boundary  is  not  perfectly  planar,  as  the  computer  simulations  assume.  Of  par¬ 
ticular  interest  are  regions  of  high  flux  intensity  in  the  experimental  image  that  correspond  to  features  identified  in 
the  simulations  of  intermode  conversion  events  (e.g.,  the  right  section  of  the  feature  marked  P  and  the  features 
marked  q  and  y).  Most  of  the  highly  concentrated  flux  of  this  type  is  due  to  intermode  conversion  between  the  two 
transverse  modes  (depicted  in  purple),  but  there  is  also  some  weakly  concentrated  flux  due  to  longitudinal  to  trans¬ 
verse  intermode  conversion  (depicted  in  green).  For  these  regions,  a  significant  fraction  of  the  energy  of  a  coherent 
longitudinal  wave  could  be  converted  into  transverse  waves.  The  success  of  this  CaF2  project  was  dependent  upon 
the  careful  calibration  of  the  BEDS  data  and  required  us  to  develop  new  local  expertise  in  the  use  of  BEDS  for  mate¬ 
rials  analysis.  To  capitalize  on  that  expertise,  we  plan  to  pursue  BEDS  studies  of  GaN  thin  films  on  a  variety  of  sub¬ 
strates. 

Work  at  Bowdoin  on  phonon  refraction  at  interfaces  involved  a  variety  of  systems  and,  in  particular,  non-cubic  ma¬ 
terials  such  as  sapphire-GaN  interfaces.  New  computer  routines  were  developed  to  simulate  both  refraction  and  re¬ 
flection  at  these  interfaces.  The  codes  are  suitably  general  to  handle  materials  of  any  symmetry  class.  As  shown  in 
Fig.  16,  reflection  imaging  allows  one  to  probe  phonon  states  without  requiring  large  amounts  of  material  for  the 
transmitted  layer.  Polarization  and  frequency  dependence  of  phonon  scattering  at  the  interface  can  be  inferred  by  the 
characteristics  of  reflection  caustics  in  the  substrate  layer  alone.  Thus,  the  use  of  reflected  phonon  imaging  to  study 
phonon  states  in  thin  films  and  low-dimensional  structures  was  of  particular  interest  to  the  MURI  collaboration. 
Msall’s  group  pursued  a  series  of  experiments  characterizing  the  qualities  of  reflected  phonon  pulses  from  a  variety 
of  interfaces,  including  clean  sapphire-liquid  helium  interfaces,  wafer-bonded  interfaces  and  epitaxially  grown  GaN 
on  sapphire  (samples  from  Morkoi;).  The  results  emphasize  the  need  for  high  quality  substrates  and  well- 
characterized  phonon  detectors.  The  expertise  gained  from  single  film  trials  was  useful  to  phonon  imaging  studies  of 
GaN  clusters  in  a  GaAs  matrix  (using  samples  provided  by  Goldman).  In  this  latter  case,  the  structures  have  charac¬ 
teristic  length  scales  comparable  to  the  dominant  phonon  wavelengths  used  to  probe  them,  which  promises  to  pro¬ 
vide  an  interesting  glimpse  of  phonon  interactions  in  very  small  structures.  Two  Bowdoin  students  worked  on  pro¬ 
jects  related  to  the  MURI  grant.  They  are,  John  Koster,  a  sophomore  who  worked  with  the  new  mask  aligner  on  the 
fabrication  of  small  scale  surface  patterns,  to  be  used  in  studying  the  effect  of  surface  patterning  on  phonon  reflec¬ 
tion,  and  Elizabeth  Robinson,  a  junior,  who  worked  on  the  computer  simulation  of  non-cubic  interfaces  and  is  con¬ 
sidering  a  senior  thesis  on  phonon  reflection  imaging  experiments.  The  Msall  group  at  Bowdoin  also  helped  charac- 


Fig.  15  -  Phonon  images  of  CaF2  with  a  grain  boundary.  Sample  thickness  =  2  mm,  scan  range  6  mm  x  6  mm.  a)  Computer 
simulation  limited  to  mode  conserving  scattering  at  interface  (FT— >  FT  blue,  ST— >  ST  red,  L— >  L  yellow);  b)  Computer 
simulation  limited  to  mode  mode  converting  scattering  (FT->ST  and  ST-»FT  purple,  L-»ST  orange,  L->FT  green),  c)  Ex¬ 
perimental  image. 
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terize  a  novel  phonon  detector,  fabricated  by  the  UM  group 
of  Kurdak.  The  device,  described  in  the  next  section,  makes 
use  of  the  cyclotron  resonance  in  an  extremely  long  Hall  bar 
to  produce  a  tunable,  narrow  band  detector.  The  measure¬ 
ment  of  the  response  of  this  device  to  phonon  signals  was 
undertaken  with  the  assistance  of  Professor  Dale  Syphers  of 
Bowdoin  College,  who  has  many  years  experience  in  Quan¬ 
tum  Hall  measurements.  The  coordination  of  his  12  Tesla 
magnet  with  optical  access  and  the  phonon  imaging  appara¬ 
tus  was  a  vital  component  of  these  measurements. 

2.1.6  PHONON  DETECTORS 

Hall  bars  as  phonon  detectors  -  Msall’s  group  observed  a 
very  strong  interaction  between  acoustic  phonons  and  a  2D 
electron  gas  confined  in  an  AlAs  quantum  well.  The  strength 
of  this  interaction  gives  these  quantum  wells  the  potential  to 
be  a  very  effective  phonon  detector.  In  these  experiments,  a 
small  AlAs  QW  structure  on  one  side  of  a  GaAs  substrate  is 
used  to  detect  phonons  generated  in  a  laser-heated  spot  on 
the  other  side.  The  AlAs  structure  is  shaped  like  a  cross 
measuring  30  pm  across  and  has  4  contacts  that  facilitate 
both  phonon  experiments  and  electric  characterisation. 

Figure  17(a)  shows  an  image  of  the  phonon  induced  voltage 
changes  measured  between  two  opposite  contacts  while 
scanning  the  laser  spot  over  an  area  of  about  1 .5  x  1 .5  mm2. 

Bright  and  dark  areas  correspond  to  positive  and  negative 
voltages  measured  across  the  AlAs  device,  respectively.  The 
main  features  of  the  image  resemble  those  of  the  well-known 
phonon  images  of  an  (OOl)-oriented  GaAs  crystal.  The  pho¬ 
non  signal  is  positive  and  negative  depending  on  the  direction  of  the  wave  vector  of  the  absorbed  phonons.  The  size 
of  the  phonon  signals  from  this  AlAs  structure  is  at  least  a  factor  of  10  larger  than  the  signals  in  comparable  GaAs 
structures  [13], 

Fig.  17  (b)  and  (c)  show  the  expected  phonon  signal  based  on  simulations  that  include  phonon  coupling  by  the  pie¬ 
zoelectric  interactions  and  the  deformation  potential  coupling  of  phonons  to  the  two  X-electron  valleys,  respectively. 
The  experimental  results  agree  well  with  the  deformation  potential  simulation,  particularly  the  coupling  to  the  slow 
transverse  (ST)  modes.  The  symmetry  of  the  experimental  image  features  also  indicates  that  the  valley  splitting  is 
negligible  on  the  energy  scale  of  the  phonons  (~  1  meV)  of  this  experiment.  The  experimental  image  also  shows 
strong  phonon  intensities  near  the  center  that  do  not  arise  from  either  the  deformation  potential  or  the  piezoelectric 
interaction.  These  features  may  be  a  consequence  of  the  detector  geometry.  Additional  simulations  that  include  this 
factor  are  planned.  The  same  AlAs  structures  show  very  pronounced  fractional  quantum  Hall  effects  if  they  are 
cooled  down  to  temperatures  below  100  mK.  This  makes  these  quantum  wells  an  attractive  alternative  to  GaAs  het¬ 
erostructures  for  the  phonon  studies  of  the  fractional  quantum  Hall  effect. 

Low  dimensional  electron  systems  as  phonon  detectors  -  Work  in  the  group  of  Kurdak  focused  on  the  development 
of  tunable  phonon  detectors  based  on  low-dimensional  electron  systems.  In  particular,  they  have  explored  using 
quantum  wire  and  quantum  dot  devices  as  phonon  detectors.  To  characterize  the  quantum  wire  detectors,  Kurdak’s 
group  deposited  A1  on  the  backside  of  the  GaAs  substrate.  Phonons  were  generated  by  hitting  the  A1  with  a  short  la¬ 
ser  pulse.  The  phonons  travel  ballistically  in  GaAs  at  low  temperatures.  They  were  able  to  detect  such  ballistic  pho¬ 
nons  by  operating  the  quantum  wire  detector  with  a  cryogenic  HEMT.  Alternatively,  on  some  of  the  devices,  they 
fabricated  electrical  heaters  on  the  back  side  of  the  GaAs  substrate.  We  have  been  comparing  phonon  response  of 
quantum  wire  detectors  for  the  two  cases  when  the  phonons  are  generated  electrically  and  optically.  We  find  that 
there  are  a  broad  range  of  time  scales  involved  with  the  relaxation  of  electrons  excited  by  phonons  in  such  sensors. 
The  phonons  generated  by  the  electrical  pulses  are  typically  detected  using  lock-in  techniques.  A  typical  phonon  re¬ 
sponse  measurement  of  the  quantum  wire  device  is  shown  in  Figure  18. 
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Fig.  16  -  Simulated  phonon  images  of  GaN  on  sapphire 
(substrate  and  film  of  equal  thickness).  The  images  show 
the  spatial  distribution  of  phonons  that  have  undergone  a 
particular  type  of  mode-conversion  at  the  interface 
(modes  are  sorted  by  velocity  and  labeled  L,  ST  and  FT 
for  longitudinal,  slow  transverse,  and  fast  transverse).  The 
probability  of  mode  conversion  at  the  interface  is  deter¬ 
mined  by  requiring  continuity  of  displacement  and  strain 
at  the  boundary.  In  this  figure,  all  mode  conserving  and 
mode  converting  possibilities  are  given  equal  weight. 
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Fig.  17  -  (a)  Phonon  drag  image  of  an  AlAs  quantum  well.  The  center  of  the  image  corresponds  to  the  [001]-  direction.  The  di¬ 
agonals  of  both  images  run  along  the  [100]  and  [010]  planes.  The  voltage  contacts  to  the  2DEG  device  are  at  the  top  and  the  bot¬ 
tom,  respectively,  b)  Simulation  based  on  piezoelectric  electron-phonon  coupling,  c)  Simulation  based  on  deformation  potential 
electron-phonon  coupling. 

Quantum  dot  devices  have  also  been  fabricated.  Potentially,  in  addition  to  being  tunable,  a  quantum  dot  detector  can 
have  enough  sensitivity  to  detect  a  single  phonon.  To  achieve  single  phonon  detection  limit,  quantum  dot  detectors 
must  be  operated  in  a  dilution  refrigerator.  A  SEM  picture  of  a  quantum  dot  detector  is  shown  in  Fig.  19.  To  perform 
phonon  experiments  at  dilution  refrigerator  temperatures,  we  must  generate  phonons  electrically.  Kurdak’s  groups 
recent  theoretical  calculations  indicate  that  the  phonon  response  in  such  devices  are  expected  to  much  larger  in  hole 
systems.  They  are  currently  fabricating  quantum  wire  and  quantum  dot  sensors  using  a  high  quality  GaAs/AlGaAs 
heterostructure  with  a  two-dimensional  hole  gas. 

Kurdak’s  group  has  also  been  exploring  the  use  of  double  quantum  well  structures  for  detecting  THz  LA  phonons. 
Unlike  conventional  tunnel  devices,  the  tunnel  coupling  between  two  2DEGs  can  be  very  strong  when  the  carrier 
densities  of  the  layers  are  identical.  This  particular  feature  has  been  used  to  realize  a  novel  double  quantum  well 
field  effect  transistor  by  a  group  from  Sandia  National  Laboratories.  A  similar  structure  has  also  been  used  to  detect 
600  GHz  radiation.  Collaborations  with  the  Sandia  group  have  begun  to  investigate  the  use  of  the  same  structures  as 
phonon  detectors  in  the  same  frequency  range.  Unlike  quantum  wire  and  quantum  dot  devices  such  phonon  sensors 
will  operate  at  zero  magnetic  field.  Phonon  energy  of  the  sensor  will  be  achieved  by  controlling  the  relative  Fermi 
energy  of  the  parallel  two-dimensional  electron  gases  using  a  back  gate.  They  have  fabricated  such  sensor  and 
shown  that  demonstrated  that  the  Fermi  energies  can  indeed  be  tuned  in  the  desired  range. 
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Fig.  18  -  Ballistic  phonon  response  of  a  quantum  wire  detector. 


Fig.  19  -  SEM  picture  of  a  quantum  dot  detector 
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2.2  Phonon  Engineering  and  Active  Control 

2.2.1  NONLINEAR  PROPAGATION  OF  THz  SOUND 

Motivated  by  the  experiments  conducted  by  the  group  of  Ted 
Norris  (University  of  Michigan),  Khurgin’s  group  have  pursued 
modeling  the  nonlinear  propagation  of  acoustic  pulses.  The  ex¬ 
periments  involved  laser  excitation  of  picosecond  acoustic 
pulses  and  their  propagation  in  bulk  Si.  Previously,  some  soli- 
ton-like  features  had  been  observed  and  explained.  In  the  final 
year  of  this  grant,  Khurgin’s  efforts  concentrated  on  adding  dif¬ 
fraction  to  the  picture  -  they  rigorously  derived  propagation 
equations  that  include  effects  due  to  nonlinearity,  dispersion 
and  diffraction.  For  the  first  time  the  anisotropic  diffraction 
equation  was  obtained.  Using  these  equations,  Khurgin’s  group 
successfully  modeled  the  pulse  propagation  numerically  and 
achieved  results  that  are  very  close  to  the  experimental  ones; 
see  Fig.  20  [14], 


DATA  SIMULATIONS 


TIME  DELAY  (ps) 

Fig.  20  -  Experimental  results  and  simulations 
for  picosecond  acoustic  pulses  experiencing  ef¬ 
fects  of  dispersion,  diffraction  and  nonlinearity. 


2.2.2  OPTICAL  PHONONS  IN  STRUCTURES  WITH  PERIODICALLY  INVERTED  POLAR  DOMAINS 

Khurgin’s  group  completed  their  work  on  periodically  inverted  polar  domains  -the  quasi-phase-matched  structures 
for  polar  optical  phonons  (Fig.  21).  They  calculated  the  dispersion  relations  in  these  structures  and  derived  the  inter¬ 
action  Hamiltonians  and  selection  rules  for  various  electronic  and  optical  processes  involving  phonons  in  PIPDS  - 
Frohlich  and  deformation  potential  scattering,  Raman  and  infrared.  The  results  shown  in  Fig.  22  demonstrate  that  the 
selection  rules  are  drastically  different  in  PIDPS.  These  results  are  reported  in  [15]. 
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Fig.  21  -  Periodically-inverted  polar  domain  structure. 


Fig.  22  -  Location  of  scattered  optical  phonons  in 
the  Brillouin  zone  of  a  PIPDS  structure.  The 
strength  of  each  process  is  represented  by  the 
gray  scale. 


2.2.3.  STIMULATED  POLARITON  SCATTERING  IN  INHOMOGENEOUSLY  BROADENED  MEDIA 

Khurgin’s  group  pursued  work  on  the  new  and  exciting  topic  of  stimulated  intersubband  Raman  scattering  by  polari- 
tons.  Fig.  23  (a)  compares  a  typical  three-level  intersubband  laser  with  the  Raman  intersubband  process;  see  Fig. 
23  (b).  Whereas  they  look  similar,  the  re-lasing  should  always  have  higher  gain  than  that  of  a  virtual  Raman  process. 
However,  experiments  by  the  NRC  group  show  that  Raman  processes  dominate  -  a  fact  that  is  difficult  to  explain. 
Khurgin  and  co-workers  developed  a  theory  in  which  the  Raman  process  does  not  end  up  in  a  given  one  -particle  in- 
tersubabnd  excitation  but  results  in  stimulated  emission  of  quasi-particles  consisting  of  intersuband  plasmon,  pho¬ 
non  and  photon.  They  show  that  the  gain  for  such  process  in  inhomogeneously  broadened  media  is  stronger  than 
lasing  gain  due  to  the  effect  of  motional  broadening. 
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Fig.  23  -  Intersubband  Processes:  (a)  lasing,  (b)  stimulated  Raman  and  (c)  stimulated  Raman  polariton 


2.2.4  OPTIMIZED  PHONON  ENGINEERING  OF 
NONLINEAR  OPTICAL  DEVICES 

The  group  led  by  Khurgin  worked  on  the  design  of  mul¬ 
tiple  GaN/AIN  QW  structures  for  nonlinear  optical  de¬ 
vices  at  communication  wavelength  (1550  nm).  The 
purpose  of  this  study  was  to  develop  and  optimize  the 
minimum  switching  power  P^-halax  where  a  is 
the  absorption  cross-section  and  i  is  the  relaxation  time 
at  a  given  bandwidth  B.  Intersubband  transitions  in  ni¬ 
tride  QWs  have  large  absorption  cross-sections  (>  1015 
cm'2)  but  the  relaxation  times  are  at  least  an  order  of 
magnitude  smaller  than  necessary  for,  say,  40GBs 
bandwidth.  The  relaxation  times  are  determined  by  LO 
phonon  and  piezoelectric  scattering.  They  showed  that 
by  using  coupled  QW’s  (Fig.  24)  one  can  tailor  the  re¬ 
laxation  time  to  the  bandwidth  to  get  the  best  perform¬ 
ance. 

2.2.5  QUANTUM  DOT  LASERS 

Understanding  of  carrier  dynamics  and  hot-carrier  effects  in  quantum  dots  and  novel  design  of  the  carrier  injection 
scheme  has  led  to  self-organized  QD  lasers  with  performance  characteristics  comparable  to  or  better  than  those  of 
QW  lasers.  Benefiting  from  the  large  density  of  states,  quantum  dot  lasers  demonstrate  large  T0,  low  threshold  cur¬ 
rent,  large  modulation  bandwidth  and  low  chirp.  For  open-air  applications,  there  is  a  need  for  reliable  lasers  which 
are  also  eye-safe.  This  necessitates  operation  wavelengths  in  the  range  of  1.3-1. 55pm,  which  is  also  the  range  of 
preference  in  fiber  optical  communication  systems.  While  great  progress  has  been  made  with  devices  emitting  at 
~lpm,  much  less  has  been  achieved  with  devices  emitting  at  1.3  and  1.55pm.  Furthermore,  there  are  unique  prob¬ 
lems  that  limit  the  performance  of  conventional  SCH  QD  lasers,  compared  to  what  is  expected  from  “ideal”  QD  la¬ 
sers  with  near  singular  density  of  states.  Bhattacharya’s  group  proposed  the  development  of  high  performance  quan¬ 
tum  dot  lasers  emitting  at  1.3  pm.  Specifically,  two  solutions  were  proposed  to  circumvent  the  inherent  unique  prob¬ 
lems  of  conventional  SCH  QD  lasers,  namely,  tunneling  injection  (TI)  and  p-doping  of  the  dots.  In  the  following,  we 
summarize  the  achievements  and  progresses  in  1.3pm  p-doped  and  tunnel  injection  lasers. 

1.3pm  p-Doped  Quantum  Dot  Lasers  -  The  p-doped  SCH  QD  laser  heterostructures  were  grown  by  MBE  on  (001) 
GaAs  substrate.  Both  undoped  and  p-doped  QD  lasers  were  grown  and  fabricated.  In  the  doped  lasers,  the  modula¬ 
tion  doping  of  the  dots  was  achieved  by  delta-doping  the  GaAs  barrier/waveguide  regions  separated  from  the 
neighboring  quantum  dot  layer  by  14  nm.  The  doping  concentration  was  varied  to  provide  sheet  acceptor  concentra- 
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Fig.  24  -  QW  with  phonon-engineered  optimal  relaxation  time. 
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tion  per  quantum  dot  layer  varying  in  the  range  (0-2)*  10'2  cm'2.  The  optimum  doping  level  was  determined  by 
studying  the  luminescence  of  the  dots  and  the  device  characteristics.  From  the  LI  characteristics  of  the  broad  area  la¬ 
sers  of  varying  cavity  length  at  15°C,  we  determine  the  value  of  internal  quantum  efficiency,  t]j,  and  cavity  loss,  y,  to 
be  62%  and  6.6  cm'1,  respectively,  for  the  p-doped  samples  with  an  optimized  delta-doping  of  5*10ncm'2.  The  value 
of  J,h  is  390 A/cm2  for  400pm  cavity  length.  For  the  undoped  QD  lasers,  //„  y,  and  J,h  are  0.89, 4.3  cm'1,  respectively. 
The  threshold  current  density  is  about  155  A/cm'2  for  400pm  cavity  length,  which  is  less  than  the  p-doped  samples. 
The  LI  characteristics  of  a  single  mode  p-doped  laser,  measured  at  various  temperatures,  are  shown  in  Fig.  25(left). 
Plotted  in  Fig.  25  (right)  are  the  threshold  current  and  differential  efficiency  as  a  function  of  temperature,  as  derived 
from  the  data  of  Fig.  25(left).  It  is  evident  that  /,/,  is  independent  of  temperature  in  the  range  5-75°C  ( T0=oo )  and  so  is 
the  differential  efficiency.  This  is  the  first  time  that  a  semiconductor  laser  has  displayed  such  temperature  invariant 
features.  In  contrast,  the  conventional  undoped  QD  lasers  exhibit  7V=69K  in  the  same  temperature  range.  We  attrib¬ 
ute  the  zero  temperature  dependence  of  the  p-doped  devices  to  the  significant  role  of  Auger  recombination  in  1.3  pm 
QD  lasers,  and  its  decrease  with  temperature  due  to  the  temperature  dependence  of  Auger  coefficient. 


Current  (mA)  Temperature  (°C) 

Fig.  25  -  Left:  Pulsed  light-current  characteristics  of  1.3pm  p-doped  single  mode  QD  lasers  at  different  temperatures. 
The  inset  shows  the  output  spectrum  at  15°C.  Right:  variation  of  threshold  current  and  slope  efficiency  of  the  lasers 
with  temoerature.  The  lines  are  fits  to  the  exnerimental  data. 


1.3pm  Tunnel  Injection  Quantum  Dot  Lasers  -  Bhattacharya’s  group  has  previously  demonstrated  high- 
performance  1 .0pm  TI  QD  lasers  with  ~25GHz  modulation  bandwidth,  zero  linewidth  enhancement  factor  and  neg¬ 
ligible  chirp.  This  experience  has  been  utilized  in  design,  growth  and  fabrication  of  1.3pm  TI  QD  lasers.  The  sche¬ 
matic  of  the  designed  laser  is  shown  in  Fig.  26(a).  The  active  region  consists  of  five  layers  of  capped  InAs  quantum 
dots  with  40  nm  buffers  in  between.  Recent  calculations  of  optical  phonon  assisted  transitions  in  an  asymmetric 
double  quantum  well  heterostructure  suggest  that  tunneling  can  be  more  efficient  if  the  populated  states  of  the  injec¬ 
tor  well  are  energetically  aligned  with  a  higher  sublevel  in  the  dot.  Therefore,  unlike  our  1.0pm  TI  laser  design 
where  electrons  are  injected  into  the  ground  state  of  the  dots  by  phonon-assisted  tunneling,  electrons  are  injected 
into  the  first  excited  state  from  a  9.5nm-thick  In0,27Gao  73  As  layer;  see  Fig.  26(b).  This  reduces  the  strain  of  the  epi¬ 
taxial  layers  in  the  1.3pm  heterostructure.  The  1.3  pm  tunnel  injection  QD  laser  heterostructures  were  grown  by 
MBE  on  (001)  GaAs  substrate.  Single-mode  ridge  waveguide  lasers  (3-5pm  ridge  width)  were  fabricated  by  stan¬ 
dard  lithography,  wet  and  dry  etching,  and  metallization  techniques.  200-2000pm  long  lasers  were  obtained  by 
cleaving.  The  light-current  characteristic  of  the  lasers  at  15°C  is  presented  in  Fig.  27.  The  threshold  current  is  less 
than  25mA  for  5pm*  1500pm  cavities.  The  spectrum  of  the  device  lasing  at  about  1255  nm  is  also  shown  in  Fig.  27. 
Measurement  of  modulation  bandwidth  and  other  dynamic  properties  of  the  lasers  are  currently  underway.  It  is  ex¬ 
pected  that  1.3pm  TI  QD  lasers  exhibit  >15GHz  modulation  bandwidth,  with  ultra-low  chirp  and  linewidth  en¬ 
hancement  factor,  which  will  be  a  breakthrough  in  the  performance  characteristics  of  1.3pm  light  sources  for  open- 
air  and  other  hosts  of  applications. 
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Fig.  26  -  (a)  The  heterostructure  of  1.3|im  tunnel  injection  QD  laser;  (b)  Band  diagram  of  the  active  region,  show¬ 
ing  tunneling  into  the  first  excited  states  of  the  dots. 


2.2.6  PHONON  CONTROL  FOR  ENHANCED  PERFORMANCE  OF  THz  SOURCES 

The  goal  of  this  portion  of  the  project  is  to  investigate  transport  and  injection  properties  of  semiconductor  structures 
in  order  to  realize  RF  sources  that  will  operate  at  frequencies  between  several  hundred  GHz  and  several  THz.  This 
range  corresponds  to  RF  periods  between  several  picoseconds  and  a  fraction  of  a  picosecond.  Details  of  the  transport 
at  this  time  scale  are  needed  to  optimize  device  operation.  Transport  properties  are  determined  by  electron-phonon 
interactions  that  are  best  studied  using  Monte  Carlo  techniques  instead  of  techniques  that  depend  on  averages  over  a 
larger  number  of  scattering  events.  It  is  also  important  to  investigate  both  high  mobility  materials  such  as  GaAs  and 
wide  bandgap  materials  such  as  GaN.  Finally,  we  need  to  study  the  injection  properties  of  tunneling  contacts  and  the 
details  of  ionization  on  short  time  scales  in  order  to  optimize  the  device  performance.  This  report  will  briefly  de¬ 
scribe  the  analysis  of  tunneling  injectors  for  transit  time  devices,  describe  the  modeling  of  carrier  drift  for  THz  time 
scales  and  finally  use  these  results  to  study  device  operation  at  THz  frequencies. 


Fig.  27  -  Pulsed  light-current  characteristics  of  1.3pm 
tunnel  injection  QD  lasers  at  15°C.  The  inset  shows  the 
output  spectrum. 
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Fig.  28  -  Conduction  band  for  biased  hetero¬ 
junction  barrier  tunneling  injector. 
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Tunnel  Injectors  -  Haddad’s  group  has  carried  out  a  detailed  study  of  tunneling  injection  into  semiconductors  using 
heterojunction  barriers.  The  band  configurations  of  the  proposed  structure  are  shown  in  Fig.  28.  The  zero  bias  condi¬ 
tion  is  shown  at  the  top  with  the  yellow  portion  above  the  contact  Fermi  level  when  bias  is  applied.  The  barrier 
changes  from  nearly  rectangular,  to  triangular  with  a  width  equal  to  the  physical  barrier  width,  to  a  triangular  barrier 
that  is  much  narrower  than  the  physical  barrier.  The  resulting  current  vs.  voltage  characteristics  for  a  100  A  GaAs- 
AlGaAs  barrier  are  shown  in  Fig.  29.  The  injection  properties  of  the  barriers  are  independent  of  the  barrier  thickness 
for  all  but  the  very  narrow  barriers.  The  transport  properties  also  depend  on  the  barrier  height,  which  in  this  system 
depends  on  the  barrier  A1  fraction.  Transport  properties  for  a  range  of  barrier  A1  fractions  are  shown  in  Fig.  30.  Al¬ 
though  the  A1  fraction  effects  the  turn  on  electric  field,  the  junction  conductance  at  a  given  current  density  is  ap¬ 
proximately  the  same.  The  modest  current  densities  and  conductance  of  these  tunneling  barriers  limits  the  THz  op¬ 
eration  of  tunneling  based  injectors.  These  limitations  will  be  discussed  after  first  describing  the  transient  transport. 

Monte  Carlo  Based  Investigation  of  Transient  Transport  -  The  properties  of  electron  transport  over  small  distances 
and  short  time  scales  can  be  investigated  using  Monte  Carlo  techniques.  We  have  modified  several  Monte  Carlo 
programs  obtained  with  the  help  of  Prof.  Umberto  Ravaioli  of  The  University  of  Illinois.  Typical  calculation  results 
for  GaAs  are  shown  in  Fig.  31 .  This  figure  shows  that  for  time  scales  on  the  order  of  50  fsecs  injected  electrons  are 
moving  nearly  an  order  of  magnitude  faster  than  the  steady  state  velocity.  Similar  transient  response  occurs  in  sili¬ 
con  and  GaN.  This  transient  “ballistic”  transport  can  be  used  for  pure  ballistic  transient  devices  operating  at  fre¬ 
quencies  above  several  THz.  However,  these  devices  need  to  operate  at  very  low  voltages  and  produce  limited 
power.  The  overshoot  degrades  more  conventional  TUNNETT  operation.  A  brief  summary  of  this  operation  is  given 
in  the  next  section. 


Fig.  29  -  Tunneling  properties  for  different  barrier  widths. 


THz  TUNNET  Device  Operation  -  The  injection  calculations  and  transit  transport  results  were  combined  to  study 
the  THz  performance  of  TUNNETT  devices.  Typical  waveforms  for  a  device  at  320  GHz  are  shown  in  Fig.  32.  The 
red  curve  is  the  voltage  across  the  structure,  the  green  curve  is  the  induced  current  based  on  the  Monte  Carlo  tran¬ 
sient  simulation  and  the  blue  curve  is  the  induced  current  based  on  an  equilibrium  transport  approximation.  The  goal 
in  a  transit  time  device  design  is  to  reduce  the  induced  current  in  the  first  half  of  the  RF  period  and  to  increase  it  in 
the  second  half.  The  combination  of  the  overshoot  and  the  tunneling  injection  is  having  the  opposite  effect.  Even 
though  the  transient  overshoot  time  is  a  small  fraction  of  an  RF  period  it’s  effect  is  spread  out  over  most  of  the  first 
half  of  the  RF  cycle  by  the  tunneling  barrier.  The  larger  current  between  n/2  and  n  phase  degrades  the  power  avail¬ 
able  from  the  device.  Similar  effects  occur  for  other  material  systems. 
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Ballistic  Operation  -  As  was  reported  in  [16]  and  12 
through  subsequent  calculations,  we  have  shown 
that  the  frequency  characteristics  of  the  THz-range  1 0 
ballistic  transit-time  oscillator  with  a  tunnel  elec- 
tron  emission  (ballistic  TUNNETT-oscillator)  sug-  E  g 
gested  in  the  same  works  [16,17]  does  not  only  de-  ^ 
pend  substantially  on  parameters  of  the  transit  space  ~  g 

(drift  region)  but  also  specifically  on  the  AC  com- 
ponent  of  the  tunnel  emission  current.  If  frequen-  o  ^ 
cies  are  low,  this  AC  component  is  defined  by  the  75 
differential  transparency  of  the  tunnel  emitter,  op  '> 
and  does  not  depend  on  the  frequency.  But  it  is 
known  [18,19]  that  such  a  quasistatic  approach 
may  be  unsuitable  in  the  THz  range  because  the  so-  q 

called  time  of  tunneling  becomes  comparable  to  the 
AC  period  or  even  exceeds  it.  Inadequacy  of  the  Time  (ps) 

quasistatic  approach  can  be  confirmed  by  our  recent 

article  [20]  where  we  have  developed  an  analytic  Fi6-  31 ' Transit  resP°nse  for  electrons  in  GaAs- 

theory  of  the  ballistic  TUNNETT,  which  can  be 

constructed  in  the  asymptotic  case  of  the  saturated  nonparabolic  dispersion  relation  for  electrons:  e(p)  =  Vs(p-ps ) 
where  Vs  is  a  saturated  ballistic  velocity.  In  this  case,  we  can  display  a  certain  connection  between  the  optimal  val¬ 
ues  ofas(a>)  (that  optimize  the  quality  factor  of  the  device)  and  the  frequency  co:  <rs{w) »  cokd  (a)  where  kd(o>)  is  a 
dielectric  constant.  To  provide  such  a  connection  in  the  quasistatic  case,  we  have  to  thin  our  tunnel  barrier  with  an 
increase  in  co.  This  thinning  leads  to  an  increase  in  the  direct  current  component  and  to  a  decrease  in  the  device  effi¬ 
ciency.  The  situation  can  be  changed  in  the  case  when  the  quasistatic  approach  fails  and  the  sufficiently  large  values 
of  cri.(ffl)can  be  reached  by  using  the  AC  tunneling  with  the  virtual  absorption  of  the  AC  photon  having  the  energy 
ha  .  Such  an  effect  [21  ]  takes  place  in  the  intermediate  region  of  the  frequencies  and  the  amplitudes  of  the  AC  elec¬ 
tric  field.  The  tunnel  barrier  is  too  thick  for  the  direct  tunneling  of  the  electron  with  the  Fermi  energy  e;.  but  is  much 
more  transparent  for  electrons  with  the  energy  eF  +  hco  ,  With  a  subsequent  increase  in  the  AC  field  intensity,  an  in¬ 
crease  in  the  DC  tunnel  current  with  the  real  (not  virtual)  hco  photon  absorption  leaves  behind  the  above-mentioned 
virtual  process  and  the  efficiency  falls  again. 


Future  Work  -  The  results  in  this  report  show 
the  limitations  on  the  power  available  from 
TUNNETT  devices  due  to  the  physics  of  tunnel¬ 
ing  and  transient  transport.  However  these  simu¬ 
lation  tools  can  be  used  to  overcome  these  limi¬ 
tations  and  improve  the  operation.  The  tunneling 
injection  limitation  can  be  overcome  by  design¬ 
ing  for  a  small  amount  of  avalanche  generation 
in  the  structure.  This  will  increase  the  current 
density,  narrow  the  injected  pulse  and  delay  the 
peak  phase.  All  of  these  effects  will  improve  the 
performance.  The  avalanche  will  have  little  ef¬ 
fect  on  the  transient  carrier  overshoot.  However 
proper  design  of  the  injector  should  move  the 
peak  of  the  induced  current  shown  in  Fig.  32 
from  before  to  after  the  n  phase  point.  The  same 
physical  effect  that  was  degrading  the  perform¬ 
ance  in  the  pure  tunneling  mode  will  be  improv¬ 
ing  it  for  the  mixed  mode  design.  It  can  be  seen 
from  the  simulation  results  that  phonon  scatter¬ 
ing  actually  aids  the  performance  in  such  de¬ 
vices  and  will  help  in  identifying  proper  material 
systems. 


Time  (ps) 

Fig.  32  -  RF  voltage  and  induced  current  waveforms  for  320  GHz  GaAs 
TUNNETT. 
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2.2.7  PHONON  ENGINEERING  IN  Si  AND  GaAs  HETEROSTRUCTURES:  LIGHT  EMITTERS 

During  this  program,  research  at  Wang’s  group  focused  on  Si  light 
emitting  diodes.  The  schematic  configuration  of  their  Si  LEDs  is  shown 
in  Fig.  33.  Three  photolithography  processes  were  used  to  fabricate  the 
Si  LEDs.  Their  processing  began  with  a  wet  oxidation  on  an  n-type  Si 
wafer  at  1000  °C  for  30  min  to  grow  5000  A  wet  oxide.  After  pattern¬ 
ing  the  active  areas  of  0.15  x  0.15  cm2  for  electroluminescence  meas¬ 
urements,  a  boron  implant  at  an  energy  of  30  keV  and  a  dose  of  1 .0  x 
1015  cm2  was  carried  out,  followed  by  an  anneal  at  600  and  1 100  °C  for 
60  s  in  an  RTP  610  rapid  thermal  processor  in  an  N2  ambient.  We  used 
aluminum  in  direct  contact  with  the  implanted  active  areas  as  the  p-side 

contact.  The  A1 
was  sputtered  in  a 

CVC  601  sputtering  system  at  room  temperature  with  a  vac¬ 
uum  of  2x1 0‘6  Torr  prior  to  deposition.  The  backside  electrode 
was  made  of  Au,  which  was  lift  emission  of  the  LEDs  an¬ 
nealed  at  600  and  700  °C  is  intermediate  between  those  an¬ 
nealed  at  800  and  900  °C,  and  those  at  1000  and  1100  °C.  Fig. 
34  compares  the  EL  characteristics  of  the  LEDs  annealed  at 
800  and  900  °C  at  currents  of  5,  7  and  10  mA.  For  the  three 
currents,  we  see  that  the  emission  of  LEDs  annealed  at  800  °C 
is  always  higher  than  that  of  LEDs  annealed  at  900  °C.  In 
summary,  we  fabricated,  characterized  and  investigated  Si 
LEDs.  We  found  optimum  processing  conditions  leading  to 
the  significant  emission  intensity  enhancement.  It  is  expected 
that  their  work  on  fabrication  and  characterization  of  Si  LEDs 
will  be  continued  so  that,  by  the  end  of  2006,  they  plan  to 
achieve  a  further  increase  in  the  EL  intensity  by  the  optimiza¬ 
tion  of  Si  LED  structure  and  fabrication  process.  Research  will 
be  addressed  to  the  doping  profiles  optimization,  growth  con¬ 
dition  improvement,  and  after  growth  temperature  treatment. 
They  also  plan  to  demonstrate  phonon-assisted  luminescence 
enhancement  in  Si/Ge/Si  tunnel  diodes  and  in  the  modified  Si 
Esaki  diodes  as  well. 

2.2.8  DILUTE  NITRIDE  GROWTH  FOR  PHONON  ENGINEERING 

During  the  considered  period,  Goldman’s  group  pursued  investigations  of  the  epitaxial  growth  of  GaAsN  and  In- 
GaAsN  using  their  Reactive  Molecular  Beam  Epitaxy  (RMBE)  system,  towards  the  development  of  advanced  het¬ 
erostructures  for  phonon-dispersion  engineering.  Goldman’s  group  made  significant  progress  in  optimizing  the 
growth  of  these  alloy  films  which  for  the  development  of  InGaAsN-based  heterostructures  for  phonon  dispersion 
engineering.  In  particular,  correlations  were  examined  between  the  growth,  structure,  electronic,  and  optical  proper¬ 
ties  of  InGaAsN  films  and  heterostructures  grown  by  solid-source  molecular  beam  epitaxy  using  an  N2-RF  plasma 
source. 

Epitaxial  Growth:  Nitrogen  Plasma  Chamber  -  With  partial  support  of  this  MURI  program,  Goldman’s  group  up¬ 
graded  their  MBE  system  to  include  a  separately  valved  and  pumped  chamber  for  the  rf  plasma  source.  The  avail¬ 
ability  of  this  plasma  chamber  enables  them  to  grow  heterostructures  with  controlled  nitrogen  incorporation  in  suc¬ 
cessive  layers,  a  capability  not  available  in  most  MBE  systems  equipped  for  growth  of  GalnNAs  and  related  materi¬ 
als.  Figures  34(a)  and  (b)  show  the  old  and  new  configurations  of  the  rf  plasma  source.  In  the  old  configuration, 
shown  in  Fig.  34(a),  which  is  typical  of  most  MBE  systems,  the  nitrogen  flux  is  controlled  by  a  standard  shutter;  in 
this  case,  nitrogen  incorporation  into  GaAs  films  occurs  even  with  the  shutter  closed.  For  example,  for  a  flux  which 
results  in  3%  N  in  GaAsN,  0.1%  N  is  incorporated  when  the  shutter  is  closed,  as  revealed  by  the  x-ray  rocking 
curves  (not  shown).  In  the  new  configuration,  shown  in  Fig.  34(b),  the  RF  plasma  source  is  contained  in  an  inde¬ 
pendently  turbo-pumped  chamber  attached  to  the  MBE  source  flange  via  a  pneumatic  gate  valve.  In  this  case,  the 
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Fig.  34  -  EL  data  of  Si  LEDs  (a)  EL  data  at 
constant  current  (7  mA)  (b)  EL  characteristics 
of  the  LEDs  annealed  at  800  and  900  °C  at  cur¬ 
rents  of  5,  7  and  10  mA. 


Fig.  33  -  Schematics  of  the  silicon  LED. 
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Fig.  34  -  (a)  original  plasma 
source  configuration,  with  nitro¬ 
gen  flux  controlled  by  a  shutter; 
(b)  plasma  source  configuration 
with  separately  pumped  and 
valved  plasma  chamber. 


pneumatic  gate  valve  is  used  to  control  the  nitrogen  flux,  such 
that  layers  may  be  grown  with  and  without  nitrogen.  This  al¬ 
lows  the  reseracher  to  grow  well-controlled  and  high  quality 
multilayers  and  heterostructures,  resulting  in  record  electron 
mobility  values  for  modulation-doped  heterostructures,  as  will 
be  discussed  below. 

Controlling  the  Nitrogen  Composition-  For  low  N  composi¬ 
tions  in  InGaAsN,  many  groups  have  found  that  the  group  V 
sticking  coefficients  are  close  to  unity,  suggesting  that  the  N 
incorporation  is  inversely  proportional  to  the  group  III  growth 
rate  [22],  For  example,  Fig.  35  presents  a  plot  of  N  composi¬ 
tion  vs.  the  inverse  GaAs  growth  rate,  for  a  variety  of  plasma 
source  configurations  and  source  gases.  The  different  symbols 
in  the  plot  are  distinguished  by  the  N  partial  pressure  (deter¬ 
mined  by  the  RGA  described  above).  For  both  the  original  and 
modified  plasma  source  configurations  [Figs.  34(a)  and  34(b)], 
the  data  fits  a  linear-least  squares  fit  fairly  well,  consistent 
with  the  nearly  unity  group  V  sticking  coefficients.  In  the  case 
of  the  new  source  configuration  with  N2/Ar  gas  mixtures, 
shown  as  the  blue  triangles  in  Fig.  35,  significantly  lower  nitrogen  incorporation  was  observed  for  growth  rates  and 
nitrogen  partial  pressures  similar  to  those  used  with  the  original  plasma  source  configuration.  In  the  modified  source 
configuration,  the  plasma  cell  was  shifted  away  from  the  substrate  position  by  nearly  30  cm,  as  shown  in  Fig.  31(b), 
perhaps  reducing  the  flux  of  active  nitrogen  radicals  impinging  on  the  sample.  When  pure  N2  is  instead  used  as  the 
source  gas,  N  incorporation  up  to  5%  is  observed,  as  shown  by  the  red  circles  in  Fig.  35. 

In-Situ  Characterization:  Forbidden  Window  of  Growth  -  Using  a  reflection  high  energy  diffraction  (RHEED)  sys¬ 
tem,  in  conjunction  with  ex-situ  measurements,  Goldman’s  group  has  recently  identified  a  “forbidden  window”  for 
MBE  growth  of  GaAsN,  in  the  temperature  range  between  -  460°C  and  5 1 0°C.  Growth  outside  the  “forbidden  win¬ 
dow”  is  likely  to  lead  to  optimum  properties.  The  reconstructions  both  inside  and  outside  the  “forbidden  window” 
are  mapped  onto  the  GaAs  surface  phase  diagram  in  Fig.  36.  It  is  interesting  to  note  that  the  “forbidden  window” 
tends  to  occur  in  the  vicinity  of  phase-boundaries  in  the  GaAs  surface  phase  diagram.  Fig.  37  (a)-(h)  shows  RHEED 
patterns  during  the  growth  of  GaAsi_xNx  layers,  collected  along  the  [110]  and  [110]  axes.  During  the  growth  of  the 
580°C  GaAs  buffer  layer,  RHEED  revealed  a  streaky  (2x4)  pattern  (not  shown).  For  growth  of  GaAs,_xNx  layers  in 
the  range  400-425°C,  a  streaky  (2x1)  RHEED  pattern  is  apparent,  as  shown  in  Figs.  36  (a)  and  (b)[23,24].  As  the 
growth  temperature  is  increased  from  460°C,  a  spotty  (1x1)  pattern  appears,  as  shown  in  Figs.  36  (c)  and  (d).  For 
higher  temperatures,  in  the  range  520-550°C,  the  pattern  becomes  a  slightly  spotty  (3x1)  pattern,  which  later  con¬ 
verts  to  a  streaky  (3x4)  pattern  at  580°C,  as  shown  in  Figs.  36  (e)  and  (f),  and  (g)  and  (h),  respectively.  The  (3x4) 
pattern  is  in  contrast  to  the  (3x3)  pattern  reported  in  some  earlier  work  [24,25], 


Fig.  35  -  N  composition  vs.  inverse  GaAs  growth 
rate,  for  various  plasma  configurations  and  source 
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Fig.  36  (top)  -  Surface  phase  diagram  for  (001)  GaAs.  Col¬ 
ored  dots  indicate  positions  on  diagram  for  GaAsN  surface 
reconstructions  in  Fig.  37  [26]. 

Fig.  37  (left)  -  RHEED  patterns  collected  along  the  [110] 
and  [1 1 0]  axes  during  growth  of  GaAsN  films  at  (a)  and 
(b)  400  -  425°C,  (c)  and  (d)  460°C,  (e)  and  (1)  520  -  550°C, 
and  (g)  and  (h)  580°C. 


Multi-Beam  Optical  Stress  Sensor  Measurements  -  Goldman’s  group  has  also  used  multi-beam  optical  stress  sen¬ 
sor  (MOSS)  measurements  to  monitor  the  relaxation  of  stress  during  the  growth  of  GaAsN  films.  These  measure¬ 
ments  allow  them  to  determine  both  stress  relaxation  mechanisms  and  fundamental  elastic  properties,  both  of  which 
are  needed  for  optimum  design  of  heterostructures  for  phonon  dispersion  engineering.  In  the  MOSS  measurement, 
shown  in  Fig.  38,  a  laser  beam  is  passed  through  an  etalon  in  order  to  produce  three  parallel  output  beams.  These 
beams  subsequently  pass  through  the  center  viewport  of  the  MBE  source  flange,  and  are  reflected  from  the  sample, 
back  through  the  viewport,  and  into  a  CCD  camera.  Monitoring  fractional  changes  in  the  spacing  between  these  re¬ 
flected  spots,  5d/d0,  or  the  mean  differential  spacing  (MDS),  enables  a  direct  measurement  of  changes  in  the  wafer 
radius  of  curvature.  The  radius  of  curvature  is  in  turn  related  to  the  stress  in  the  film,  a,  using  Stoney’s  Equation. 
Using  the  convention  that  decreasing  beam  spacing  is  negative,  and  tensile  stress  is  positive,  the  stress  x  thickness 

product,  cjhf,  is  related  to  the  MDS  as  follows: 


MBE 

chamber 


Fig.  38  -  Schematic  diagram  of  the  multi-beam 
optical  stress  sensor.  The  laser  beam  is  split 
into  multiple  beams  by  the  etalon,  which  are  re¬ 
flected  into  a  CCD  camera  that  monitors  the 
spacing  of  the  beams. 


f  Sd\  cos  or 


ah,  =  -  — 
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where  a  is  the  angle  of  incidence,  measured  with  respect  to  the 
sample  normal,  L  is  the  total  optical  path  length  from  the  sam¬ 
ple  surface  to  the  detector,  Ms  is  the  biaxial  modulus  of  the 
substrate,  and  hs  and  hf  are  the  substrate  and  film  thickness,  re¬ 
spectively.  Goldman’s  group  has  examined  the  effects  of  in¬ 
creasing  N  composition  on  the  elastic  constants  of  GaAsN  al¬ 
loys.  A  comparison  of  stresses  determined  using  in-situ  meas¬ 
urements  of  wafer  curvature  (MOSS)  with  those  determined 
from  x-ray  rocking  curves  (XRC)  using  a  linear  interpolation 
of  lattice  and  elastic  constants  is  shown  in  Fig.  39  (a).  In  both 
cases,  the  composition  was  independently  measured  using  nu¬ 
clear  reaction  analysis.  For  N  compositions  >  1 .3%,  the  XRC 
stresses  are  higher  than  the  MOSS  stresses,  presumably  due  to 
bowing  of  the  elastic  constants  in  GaAs,.xNx  : 

Cj,(GaAS|.INI)  =x  Cy(GaN)  +  (1-x)  C,j(GaAs)  -  btJx  (1-x) 

where  a  linear  interpolation  (i.e.  Vegard’s  Law)  assumes  bj,=0, 
and  non-zero  values  of  by  imply  bowing  of  the  elastic  proper- 
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q  j - , - p - , — , - , -  ties.  Figure  39  (b)  shows  the  composition  dependence  of 

tfCRCtr  |  it  Cn  and  C12  using  for  which  MOSS  and  XRC  stresses  are 

^vfOSS  <7  I  ¥  equal,  with  a  range  of  b|/b2  ratios.  The  composition  de- 

0.4  : - — — -  i  *  pendence  of  the  elastic  constants  is  compared  with  values 

I  •  determined  using  bulk  Brillouin  scattering  measurements, 

j  1  ri  interpreted  using  the  density  and  refractive  index  of  GaAs. 

'S-  i  ,  i  "If  Further  experimental  and  computational  work  is  needed  to 

o  ]f  i  t  j  $  fully  understand  the  origins  of  the  N-composition- 

'%  q  2 .  J 1  j  dependent  bowing  of  the  elastic  properties. 

t  *  |  In-Situ  Scanning  Tunneling  Microscopy-  Although  the 

_  j  RHEED  and  MOSS  data  provide  average  information  on 

i"  I  the  surface  structure,  Goldman’s  group  has  also  used  addi¬ 

tional  tools  in  order  to 

0 _ , _  j  _ , _  understand  and  control 

15  '  '  '  '  the  film  formation 

^11  processes  on  the  near 

atomistic  scale.  They 
use  a  of  in-situ 

1 0  scanning  tunneling  mi- 

Is  "  — -  croscopy  (STM) 

q  "snapshots"  obtained  ^ SsH 

q  after  a  series  of 

§  ,i  _ _ # _  growth/quenching  se- 

g  ^0  •  quences.  An  example 

O  _bll/b12=15  of  in-situ  s™  Studies  —  *'  1 

|  bu/b|2=2  of  sequential  layers 

3  "Brillouin  during  the  growth  of 

01 _ _ _ , _ .  I  *  . . .  GaAsN  are  shown  in 

0  0.5  1.  1.  2.0  2.5  Figs.  40  (a)-{c).  For 

N  Composition  (%)  the  500  nm  GaAs 

Fig.  39  -  (a)  Stress  measured  by  in-situ  MOSS  and  calcu-  s^own  >n  F*S- 

lated  from  ex-situ  XRC  measurements,  using  a  linear  in-  ^  (a)>  a_ *  Pen°dicity 

terpolation  of  lattice  parameters  and  elastic  constants;  The  and  [110]  elongated 

stress  differences  in  (a)  are  used  to  determine  the  compo-  holes  are  apparent.  Af- 

sition  dependence  of  the  elastic  constants,  which  are  ter  0.1  nm  GaAsN 

compared  with  Brillouin  scattering  measurements  in  (b).  growth,  shown  in  Fig. 

40  (b),  isolated  do¬ 
mains  nucleate  on  flat - -  - 

regions  between  the  elongated  holes.  In  Fig.  40  (c),  after  0.3  nm  GaAsN 
growth,  a  nearly  continuous  film  is  observed.  The  GaAsN  films  also  possesses 
a  4x  periodicity,  with  a  lattice  parameter  slightly  larger  than  that  of  GaAs. 

These  results  suggest  that  N  atoms  tend  to  occupy  interstitial  sites  at  the  earliest 
stage  of  growth.  In  the  future,  Goldman’s  group  plans  to  use  the  feedback  from 
in-situ  STM  to  identify  the  appropriate  growth  conditions  which  will  allow  the 
tuning  of  interstitial  N  incorporation,  N  cluster  formation,  and  alloy  phase  sepa¬ 
ration,  in  both  GaAsN  and  InGaAsN. 

Ex-Situ  Characterization:  Ion  Beam  Analysis-  To  determine  the  substitutional 
and/or  interstitial  incorporation  of  N  in  InGaAsN  alloys  and  heterostructures, 

Goldman’s  group  used  a  combination  of  Rutherford  Backscattering  Spectrome¬ 
try  (RBS)  and  Nuclear  Reaction  Analysis  (NRA)  in  channeling  and  non-  Fig  40 .  In_situ  scanning  Umneling 
channeling  conditions.  This  set  of  measurements  enabled  them  to  estimate  the  micr0sc0py  images  of  sequential 

fraction  of  incorporated  N  atoms  sitting  in  split  interstitial  sites,  as  shown  in  layers  of  (a)  GaAs,  (b)  0.1  nm 

Fig.  41.  To  date,  Goldman’s  group  has  successfully  used  this  approach  to  study  GaAsN,  and  (c)  0.3  nm  GaAsN.  The 

N  incorporation  in  dilute  GaAsN.  Prior  to  their  work  [27],  conflicting  results  gray-scale  ranges  displayed  are  (a) 

had  been  reported  regarding  the  mechanism  of  N  incorporation  [22,28,29]  As  1.3  nm,  (b)  1.4  nm,  and  (c)  1.4  nm. 


N  Composition  (%) 


Fig.  39  -  (a)  Stress  measured  by  in-situ  MOSS  and  calcu¬ 
lated  from  ex-situ  XRC  measurements,  using  a  linear  in¬ 
terpolation  of  lattice  parameters  and  elastic  constants;  The 
stress  differences  in  (a)  are  used  to  determine  the  compo¬ 
sition  dependence  of  the  elastic  constants,  which  are 
compared  with  Brillouin  scattering  measurements  in  (b). 


Fig.  40  -  In-situ  scanning  tunneling 
microscopy  images  of  sequential 
layers  of  (a)  GaAs,  (b)  0.1  nm 
GaAsN,  and  (c)  0.3  nm  GaAsN.  The 
gray-scale  ranges  displayed  are  (a) 
1.3  nm,  (b)  1.4  nm,  and  (c)  1.4  nm. 
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shown  in  Fig.  42,  their  studies  revealed  significant  composition-dependent  incorporation  of  N  into  non-substitutional 
sites,  presumably  as  either  N-N  or  N-As  split  interstitials.  Furthermore,  we  identified  the  (2x1)  reconstruction,  asso¬ 
ciated  with  low  growth  temperatures,  as  the  surface  structure  which  leads  to  the  highest  substitutional  N  incorpora¬ 
tion,  likely  due  to  the  high  number  of  group  V  sites  per  unit  area  available  for- N-As  surface  exchange. 


Fig.  41  -  Ball  and  stick  crystal  schematic  for  GaAsN 
with  substitutional  N,  N-As  split  interstitial,  and  N-N 
split  interstitial.  The  white,  green,  and  black  spheres  rep¬ 
resent  Ga,  As,  and  N,  respectively. 
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Fig.  42  -  Substitutional  (red)  and  interstitial  (green)  concentrations 
vs.  total  N  concentration  for  GaAsi.xNx  with  varying  x.  The  concen¬ 
tration  of  interstitial  N  increases  with  total  N  concentration.  “See 
Ref.  28;  bSee  Ref.22;  “See  Ref.  29. 


Electronic  Transport  Properties:  Dilute  Nitride  Films  -  To  date,  Goldman’s  group  has  systematically  studied  the 
effects  of  growth  conditions  on  the  electron  mobility  of  GaAsN  films,  using  the  variations  in  nitrogen  plasma  source 
configurations  discussed  earlier.  As  shown  in  Figs.  43  (a)  and  43  (b),  the  electron  mobility  of  GaAsN  films  is 
strongly  dependent  on  the  plasma  source  gas  and  GaAs  growth  rate,  respectively.  In  particular,  the  highest  electron 
mobilities  are  apparent  for  films  grown  with  pure  N2,  using  higher  growth  rates.  As  shown  in  Fig.  43  (a),  the  elec¬ 
tron  mobilities  of  our  GaAsN  films  grown  in  optimum  conditions  are  most  often  higher  than  literature  reports  to  date 
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Fig.  43  -  (a)  Electron  mobility  vs.  N  composition  for  our  GaAsN  films  (McKay),  in  comparison  with  other  reports.  The  solid 
and  dashed  lines  are  calculated  assuming  all  N  substitutes  for  As,  The  dashed  line  includes  the  resonant  scattering  contribution 
from  N  clusters.  In  addition,  the  mobilities  of  GaAsN/GaAs:Si  SLs,  shown  schematically  in  (b),  are  included  in  (a).  The  corre¬ 
sponding  energy  band  diagram  of  the  SLs  in  (b)  is  shown  in  (c). 
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Fig.  44  -  Cross-section  of  selectively  doped  GaAsN/AlGaAs  heterostructure  (left)  and  corresponding  conduction  energy 
band  diagram  (right).  A  two-dimensional  electron  gas  is  formed  in  the  triangular  quantum  well  potential  in  the  GaAsN 
laver. 


[30,31],  The  most  striking  feature  of  this  combined  data  is  the  steep  drop  electron  mobility  with  N  composition,  a 
clear  indication  that  N  related  defects  are  the  main  source  of  scattering  in  dilute  nitrides.  Yet,  the  highest  measured 
electron  mobilities  are  still  about  a  factor  two  smaller  than  the  recent  calculations  of  electron  scattering  by  Fahy  et 
al.  which  include  resonant  scattering  from  N  clusters  [32]. 

Superlattice  Structures-  In  principle,  due  to  the  modulation  doping  of  Si  in  a  doped/undoped  superlattice  (SL) 
structure,  the  electron  mobility  of  a  SL  structure  may  be  larger  than  a  corresponding  film  with  similar  composition 
[33].  Figure  43  (b)  presents  a  GaAsN/GaAs:Si  heterostructure,  with  the  corresponding  conduction  energy  band  dia¬ 
gram,  calculated  using  the  Poisson  Solver.  For  this  structure,  the  electron  mobility  was  improved  by  the  use  of  the 
gate-valve  interface  control  in  comparison  with  the  conventional  shuttering  approach. 


GaAsN/AlGaAs  heterostructures  -  A  selectively  doped 
GaAsN/AlGaAs  heterostructure  and  the  calculated  conduction 
band  energy  diagram  of  the  corresponding  structure  are  shown  in 
Fig.  44.  In  such  a  structure  a  two-dimensional  electron  gas 
(2DEG)  is  formed  in  a  triangular  confinement  potential  in 
GaAsN.  Similar  to  GaAs/AlGaAs  heterostructures,  the  carrier 
density  of  the  electron  gas  can  be  engineered  by  varying  the 
spacer  layer  thickness.  In  collaboration  with  Kurdak’s  group, 
Goldman’s  group  designed  and  grew  a  heterostructure  with  very 
low  N  composition  (<0.05%).  Remarkably,  at  low  temperatures, 
in  this  very  low  N  composition  2DEG  structure,  we  observe  inte¬ 
ger  quantum  Hall  effect  and  well-resolved  Shubnikov  de  Haas 
oscillations  at  magnetic  fields  down  to  1  Tesla,  as  shown  in 
Fig.  45.  The  extracted  the  carrier  density  and  the  electron  mobil¬ 
ity  of  the  2DEG  are  1.13xl04  cm2 /Vs  and  1.86x10"  cm'2,  re¬ 
spectively.  The  electron  mobility  of  this  2DEG  structure  is  sub¬ 
stantially  higher  than  those  recently  reported  in  the  literature 
[22,31]  presumably  due  to  our  gate- valve  control  of  the  het¬ 
erostructure  interfaces,  described  earlier. 


Fig.  45  -  Quantum  Hall  effect  and  Shubnikov-de 
Haas  data  at  4.2  K  from  a  GaAsN/AlGaAs  het¬ 
erostructure  with  a  two-dimensional  electron 
gas. 
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